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EXECUTIVE SUMMARY 


In conjunction with a similar study in the six New England 
States, 278 lakes and streams were sampled in the nine Middle 
Atlantic States (New York, New Jersey, Pennsylvania, Maryland, 
Delaware, Virginia, West Virginia, North Carolina, and 
Tennessee). Primary objectives of the study were to document the 
sensitivity of waters in the study area to acidification, and to 
test the validity of available models and classifications for 
predicting acidification of these waters. Waters were chosen to 
be relatively free of direct disturbance by human activity in 
their watersheds, and to be at as high an elevation as possible 
while maintaining representation of all geographic areas and 
bedrock and soil sensitivity classes. Lakes comprised 60% of the 
Sites and streams 40%. Historical data for the same waters were 
gathered wherever available. 


About 49% of the unpolluted, relatively undeveloped waters 
sampled in the nine states were sensitive to acidification on the 
basis of having alkalinity less than 200 yueq/l or of having 
underlying bedrock of low acid neutralizing capacity. Alkalinity 
is probably the most useful and relatively accurate predictor of 
sensitivity to acidification. Calcite saturation index (CSI) was 
highly correlated with alkalinity, and appeared to have no 
Significant advantage as a predictor. Of the waters sampled, 32% 
were susceptible to acidification on the basis of having CSI 
greater than 3. Neither bedrock nor soil class is a particularly 
Strong predictor of sensitivity to acidification, but bedrock 
class is somewhat correlated with alkalinity of Middle Atlantic 
States waters, while soil type is not, at least when determined 
from rock and soil maps presently available. 


Comparison with available historical data indicated that 
many of the sampled waters have decreased in alkalinity in recent 
years. A similar decrease in pH was not noted, but this is to be 
expected so long as the acid neutralizing capacity (as 
represented by alkalinity) of most waters has not been completely 
depleted. However, the accuracy of the historical data could not 
be adequately evaluated to enable confident conclusions about 
temporal trends in acidification. 


Few of tne waters showed high concentrations of aluminum, 
and only one stream had aluminum above 200 yeq/l at a pH below 
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5.5, which is believed to be the range of toxicity to fish. 
However, we noted a doubling of average aluminum concentrations 
in waters with pH below 5, compared to waters with higher pH. 


Suifate levels (mean 227 yeq/l) found in the sampled waters 
were higher than those reported from some other areas where 
precipitation is acidic, and much higher than those from areas 
with nonacidic precipitation. Sulfate was not correlated with 
acidity, perhaps due to the relatively uniform deposition of 
sulfate by precipitation in the study area. 


DH was correlated strongly with alkalinity, conductivity, 
and (negative) calcite saturation index, and weakly with calcium, 
magnesium, and disturbance. Alkalinity was strongly correlated 
with calcium, pH, conductivity, and magnesium, and weakly with 


sulfate, chloride, disturbance, rain pH, and (negative) 
elevation. Neither pH nor alkalinity was strongly correlated 
with aluminum, manganese, color, total organic carbon (as 


predicted from color), or size of waters. 


Differences in water chemistry between lakes and streams 
were small, with streams in general being less acidic. Higher 
order streams had higher conductivity and pH, as expected due to 
their larger, more varied watersheds. Waters most affected by 
human disturbance had higher pH, conductivity, and chloride. 


The calcium-pH model of Henriksen predicted the pH of waters 
in the study area reasonably well. It indicated that about 28% 
of the waters were acidified. Other published models and 
Classifications (Thompson Cation Denudation Model, Altshuller and 
McBean, Sulfate-pH) that were tried predicted water pdH poorly. 


The study results indicate that a large number of relatively 
undisturbed waters in the Middle Atlantic States are either 
already acidified or are susceptible to acidification. Thus, 
parts of the Middle Atlantic States should be included on maps 
showing regions sensitive to acid precipitation. A large portion 
of the waters we sampled in the region have decreased in pH, 
alkalinity, or both in recent years. 
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INL RCDUCTION 


The precipitation that falls on the northeastern United 
tates is highly acidic. Its pH averages about 4.2 but is often 
much lower (Likens et al. 1979; National Atmospheric Deposition 
Program 1984). Research has indicated that acidic precipitation 
has had profound effects on aquatic eccsystems, particularly 
fish, in Scandinavia (Leivestad et al. 1976), Ontario, Canada 
(Beamish 1976), and the Adirondack Mountains of New York 
(Schofield 1977). Evidence of serious effects on both aquatic 
and terrestrial ecosystems in the Northeast and other parts of 
the United States is accumulating rapidly (see Haines 1981 and 
Glass et al. 1982 for reviews). 


Tne impact of acidic precipitation on aquatic ecosystems 
depends on a variety of climatic, ageologic, topographic, 
morphometric, biotic, and anthropogenic factors. Even in areas 
where effects of acidification on lakes have been pronounced (La 
Cloche Lakes, Ontario, and Adirondack Mountains, New York) there 
are many lakes chat seem unaffected, interspersed with acidified 
lakes. This phenomenon illustrates the need for indices of 
vulnerablity that can be applied to individual lakes and streams 
and that are sensitive to local conditions of rock, soil, and 
hydrology. 


If regulations or legislation controlling emissions of 
Sulfur and/or nitrogen oxides are to be sought as a means of 
reducing the adverse effects of acidic precipitation on 
ecosystems, detailed information will be required to substantiéte 
the need for such reductions. Such information would need to 
include projections of the quantity of aquatic habitat that is 
Susceptible to acidic precipitation inputs, and dose-response 
estimates for the effects of such inputs. Dose-response 
estimates will require detailed, long-term study in a variety of 
habitats subject to varying inputs of acidic precipitation. 


The present study is an attempt to improve available 
estimates of the extent of vulnerable habitat by documenting the 
sensitivity of waters in the study area to acidification, and to 
determine which, if any, of the available predictors of 


vulnerability is most appropriate for use on these waters. The 
terms "sensitivity" and "vulnerability" as used in this report 
(and most acid rain literature) refer to no or low acid 


neutralizing capacity, thus the quality of being easily 


7 
+ 


acidified. Existing surface water chemistry data were not 
Satisfactory for validating these predictors. Existing data were 
predominantly collected in high order watersheds whe 


e soils are 
thick and geological comp] f buffering 


le 
Materials likely, while acid D 
= 3 


recivicatior likely to 
affect low order, peace akes anc streams. Also, many 
existing water chemistry date Sites were selected to determine 
the effects se human activities on water quality. Human 


activity makes natural chemistry difficult tc ascertain, and 
there is almost always significant alkalinity contained in 
Gischarages from human activities (with the exception of acid mine 
drainage and a few other wastes). In gen eral, human activity 


is 
concentrated in areas that were originally selected (for 
agricultural productivity, and are thus likely to be well 
buffered. The most pervasive anthropogenic impact on alkalinity 
may de from ages cul aeal liming [JanRSOn 1984). Finally, many 
existing water hemistry data were coliected with methodology 
(e.g., colorim tri DH, fixed endpoint alkalinity titration) that 
nrodcees inaccurate results when used in low ionic st rength 
waters. us the present study was planned to collect vaiid 
present- day data on waters as unmodified as possible, located on 
low order watersheds, and to use best available methods re low 
ionic strength waters. The data collected, in additio ae) 
indicating current degree cf acidification, if any, were ae to 
evaluate several proposed indices of vulnerability. 


The study was divided into two parts. The six New England 
States (Maine, New Hampshire, Ver Massachusetts, 
Connecticut, and Rhode island) were sam ol qd and the results 
prepared by Dr. Terry A. Haines and associates at the U. S. Fish 
and Wildlife Service field Station at the University of Maine in 
Orono, Maine (Haines and Akielaszek 1983: hereafter referred to 
as "the New England study"). We sampled sites in the nine Middle 

tlantic States (New York, New Jersey, Pennsylvania, Maryland, 
Delaware, Virginia, West Virginia, and tne Appalachian parts of 
North Carolina and Tennessee). 

Virtually identical field and laboratory methods were used 
by the two groups, and statistical analyses of the two sets of 
data were similar. In general, the two reports follo a common 


outline. 











MATERIALS AND METHODS 


The following description refers to the Middle Atlantic part 
of the study but is generally true of the New England part as 
well. Methods differing significantly between the study parts 
are footnoted. 


SAMPLING METHODS 


Selection of Sampling Sites 





Sites included in this study were selected from all 
available surface waters with respect to several criteria: 


- broad geographic coverage 
- bedrock acid neutralizing capacity 
- soil cation exchange capacity 
- location in upper portion of watershed 
(i.e., altitude, stream order) 
- absence of direct human disturbance in watershed 
- availability of historical data 
- difficulty of access 


The waters sampled were primarily lakes (168 of 278; 60%) 
(for both studies, 362 of 504 or 72%). Streams were usually 
Sampled when a suitable lake could not be located in a desired 
area, Or was not accessible. In general, the lowest order stream 
feasible was chosen. Similarly, impoundments (defined as any 
water body whose level is controlled or fixed by a man-made 
structure) were sampled only when a natural lake was not 
available in the same area. This was a somewhat arbitrary 
distinction in some cases, since many impoundments where the dam 
had no provision for water level adjustment and the level had 
apparently been stable for several years were equivalent to 
"natural" lakes. Impoundments tended to be more predominant in 
the southern part of the study area, probably reflecting the lack 
of glaciation there. Impoundments which showed evidence of water 
level fluctuation or adjustment were avoided. Lakes sampled were 
Classified as either "seepage" lakes with no visible inlet (n=17) 
or “drainage" lakes with well-defined inlets (n=151). 
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Sampling sites over the region 


Geographic distribution o t 
(Figure 1), Dut a few areas were 


was as uniform as possible 
inadequately represented. Or example, there are few undisturbed 
lakes or streams in central New York. 
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Sampling sites were selected to represent ail classes of 
bedrock acid neutralizing capacity and soil cation exchange 
Capacity present in the region. Bedrock acid neutralization 
Capacity classes were those described by Hendrey et al. (1980), 
and soil cation exchange capacity classes were those of McFee 
(1980) (Table l). 

The original maps depicting the distribution of these 


Classes were obtained from the authors (Stephen Norton, 
University of Maine, Orono, and William McFee, Purdue University, 
Lafayette, IN, personal communications). These maps were 


generally of 1:500,000 scale, the same as the U. S. Geological 
Survey (USGS) state base map. 


Waters were selected to represent the various bedrock and 
SOil classes in approximate proportion to the area of the classes 
in each state. Bedrock classes were generally given preference 
in cases where a choice nad to be made. Some waters may have 
been misclassified as to bedrock or soil class due to difficulty 
in precisely locating the water on the state base map (many small 
waters are not shown on these maps, and locations were 
transferred from 7.5 or 15 minute USGS quadrangle maps). Also, 
“smoothing” and other removal of complexity were performed in the 
creation of the bedrock and soil maps. In general, sites were 
selected at the highest elevations possible within chosen 
bedrock/soil type areas, since higher elevation, lower order 
waters are usually more susceptible to acidification and less 
likely to be affected by other forms of pollution. Héeadwater 
Sites were generally low in direct human disturbance. 


For each state, appropriate agencies were asked to supply 
existing historical water chemistry data on the chosen waters. 
Particular attention was given to selecting sampling sites for 
which there were historical data readily available; nearly half 
(123 of 278, 44%) of the final sites had such data (for both 
Studies, 218 of 504 or 43%). 


Finally, accessibility of the sites was considered in making 
a final selection. Preferred sites were at least one-half mile 
from the nearest vehicle access point, but not more than 2 miles 
from that point. Most sites had at least a marked foot trail 
leading from the access point to the site. 
















Site Type 


Lake /\ 

















Stream 





Impoundment <> 




















Figure 1. Map of the Middle Atlantic States showing sampling 
Sites and type of water body for each. 
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Table 1. Bedrock and Soil Sensitivity (to acidification) 
Classes (See Hendrey et al. 1980 and McFee 1980 
for further details) 











Code Class Description 
(Rock ) 
1 I No to low acid neutralizing capacity 
2 II Low to medium acid neutralizing capacity 
3 III Medium to high acid neutralizing capacity 
4 IV "Infinite" acid neutralizing capacity 
5 V Unknown but probably moderate to high 


acid neutralizing capacity 


(Soil) 
l Sl Sensitive soils dominate the area 
2 S2 Sensitive soils are significant, but 


cover less than 50% of the area 


3 SSl Slightly sensitive soils dominate 
the area 
4 SS2 Slightly sensitive soils are 


Significant, but cover less 
than 50% of the area. 


5 NS The area contains mostly nonsensitive 
soils 
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Sample and Field Data Collection 





Due to the large size of the study area, the sampling period 
extended over an entire year, July 1980 to August 1981, with most 
samples collected in spring and summer 1981. Some values for a 
particular water may thus reflect seasonal limnological 
conditions. 


Standard field procedures were similar to those used in New 
England. Most water bodies selected for survey were visited by a 
two person survey team. All were reached by driving as close as 
possible with a four-wheel drive vehicle, then carrying equipment 
the remaining distance. The major vegetation, exposed bedrock, 
inflow and outflow streams, aquatic vegetation, and evidence of 
human activity (roads, campsites, seasonal dwellings, permanent 
dwellings, logging, agriculture) visible from the water and 
observed on the trip into and out of the area were noted. 
Although waters with human disturbance in the watershed were 
avoided, a higher proportion of slightly disturbed waters had to 
be included in the Middle Atlantic States than in New England due 
to the scarcity of completely undisturbed sites. An arbitrary 
code was devised to represent differing levels of disturbance; 
waters were assigned a code number based on a visual survey of 
the shoreline and watershed (Table 2). A canoe was used to reach 
the deepest area of lakes and ponds except in a few cases when 
thin ice prevented access by canoe or foot beyond = shore 
Structures (@.g., ends of docks). In such cases-~ only 
near-surface samples were collected, as far from shore as 
feasible. Adequately frozen lakes were sampled by spudding a 
hole through the ice near the deepest area. "Surface" samples in 
these cases were collected by first gently removing the water-ice 
mixture from the hole and allowing it to refill, then placing the 
bottle about 0.5 m below the water surface. Water samples were 
collected in acid-washed, distilled-water-rinsed linear 
polyethylene bottles. 


In streams, water samples were collected from mid-depth and 
midstream directly into the bottles. In lakes, one set of 
samples was collected at 0.1 m below the surface directly into 
the bottles and one set near the bottom (without disturbing 
sediment) in a plastic Van Dorn-type water sampier. Temperatures 
at sample depths were measured directly with a thermistor, and a 
Secchi disc reading was recorded in lakes. Width of streams was 
estimated visually at the site. Three 500 ml and two 125 ml 
bottles of sample were collected at each depth. All bottles were 
rinsed three times with the water to be sampled before filling. 
Bottles were filled to the brim to exclude air, thus minimizing 
Gas exchange and consequent chemical changes. Immediately on 
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Table 2. Criteria Used to Estimate Relative Amount of 
Human Disturbance 





Code bisturbance 





1 No visible disturbance: access by foot trail 
Or aircraft only 


2 Slight disturbance: access by 4-wheel 
drive vehicle; one or two seasonal 
dwellings in the watershed; evidence of 
logging activity but not presently active 


3 Moderate disturbance: access by 2-wheel 
drive vehicle; more than two seasonal 
dwellings; evidence of active logging 
but not immediately adjacent to water 


4 Severe disturbance: access by paved road; 
permanent homes around lake; agriculture 
in watershed; active logging adjacent to 
water 





return to shore, 2.5 ml of 20% nitric acid was added to each 125 
ml sample to preserve it for cation analysis, and all samples 
were placed on ice. 


ANALYTICAL METHODS 


Field Analyses 





Initial measurements of pH and conductivity were performed 
in the field from 1-5 hours after collection. One 500 ml bottle 
from each sample depth was warmed to near 20 °C while still 
closed. An Orion model 399A meter equipped with an Orion model 
91-62-00 glass combination electrode (specially designed for low 
ionic. strength solutions) was used for pH determinations in the 
field*. Before each determination, the meter was standardized 
against pH 4 and 7 buffers. After being rinsed with distilled 
water and sample water, the electrode was inserted directly into 
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the bottie and moved gently until the reading of the meter 
Stabilized (about 5 minutes). This procedure minimized contact 
between sample and atmosphere, a potential source of error with 
poorly buffered waters. The procedure was repeated and both 
readings were recorded. 


Conductivity or the sample was measured directly in the 
bottle using a Presto-Tek model DP-32 conductivity/temperature 
meter. Two successive readings were recorded along with the 
associated sample temperature; the meter had a built-in 
temperature compensator which standardized readings to 25 °C. 


Laboratory Procedures 





The remaining water samples were kept on ice, returned to 
the laboratory, and frozen until analyzed~. The acidified 125 ml 
bottles were given to a subcontractor (Dr. William J. Myers, 
Micro-Technic Systems Laboratory, Pleasant Gap, PA 16823) for 
analysis of aluminum, calcium, manganese, sodium, magnesium, and 
potassium. Total aluminum was measured by the colorimetric 
method of Dougan and Wilson (1974); the other elements were 
measured by the standard methods recommended by .the JU. S. 
Environmental Protection Agency (Kopp and McKee 1979)”. 


Samples for metal analysis were not filtered, in order to 
maintain compatibility with the New England study (Haines and 
Akielaszek 1983). Some workers believe that unfiltered samples 
may be preferable for aluminum analyses of poorly buffered waters 
low in dissolved sclids (Charles Driscoll, Syracuse University, 
personal communication). All analyses were duplicated, using one 
sample from each bottle. 


2 Haines and  Akielaszek (1983) analyzed = aikalinity, pH, 
conductivity, and color in the field, and used a plastic~-body, 
gel-filled combination electrode on their portable pH meter. 


The New England samples were refrigerated, not frozen. 


Haines and Akielaszek (1983) analyzed sodium and potassium by 
air-acetylene flame atomic absorption spectrophotometry, 
calcium and magnesium by nitrous oxide-acetylene flame AAS, and 
manganese by graphite furnace AAS. They measured aluminum 
initially by the colorimetric method of Dougan and Wilson 
(1974) and later by graphite furnace AAS. The two methods gave 
“comparable results". 








The remaining two 500 mil obcttles were brought to rcom 
temperature and immediately analyzed for alkalinity, sulfate, 
chloride, and color. As a check on the field measurements, pH 
and conductivity of the thawed sample were recorded. 
Conductivity was re-measured using the same instrument and 
technique as in the field. pH was re-measured using a Corning 
mocel 125 meter equipped with the same electrode used in the 
field, and the same technique was employed. (As fewer problems 
and greater precision were noted with the lab pH measurements, 
they were used for the data comparisons which follow). 


Alkalinity was determined by the double endpoint, 
potentiometric titration method as described in American Public 
Health Association et al. (1981). Duplicate 100 ml samples were 
titrated with sulfuric acid (standardized to 0.0200 N) to pH 4.5 
and then to pH 4.2, recording the amount of acid used to reach 
each endpoint. This method has proven to be nearly as accurate 
as the Gran plot method (Stumm and Morgan 1970) and considerably 
faster. It is also more accurate than the fixed (single) 
endpoint method often used in the past by many investigators 
(Kramer and Tessier 1982). However, the values obtained were 
also used to calculate the fixed endpoint alkalinity for use in 
historical data comparisons. 


To verify the comparability of the double endpoint method 
with the Gran plot method, and the analysis of frozen samples for 
alkalinity with unfrozen samples, we collected a set of 
replicated samples from 16 waters in central Pennsylvania. Nine 
of these waters also appeared in the project samples; the 
remainder were chosen to represent a broader range of trophic and 
chemical conditions. From each water, four samples were 
collected just below the surface. Two of the four were frozen 
for 2 weeks and thawed just before analysis (similar to the 
samples of the main survey); the other two were refrigerated for 
2 days before analysis but not frozen. Samples were analyzed for 
alkalinity by both the double endpoint and Gran plot methods, and 
for conductivity, pH, sulfate, and color. A paired-comparison 
t-test (SAS Institute 1982) was used to test for differences 
between means of frozen and unfrozen samples for each analysis, 
and between double endpoint and Gran plot methods for alkalinity. 
Differences were standardized by means of the formula: 


SEGCOCONGS So S4eSessaeeesens dasa ene anaemia 
(value before + value after) / 2 
Replicates were treated as separate samples. Results are 


presented in Table 3, and show relatively small differences 
between the procedures under consideration. 
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Table 3. Comparison of Frozen and Unfrozen Sample Analyses 
and of Alkalinity Methods for 16 Waters (values 
for conductivity in uS/cm, for color in Pt-Co 
units and others in yeq/1) 











Unfrozen Mean Standardized 
Variable = seer een----- Difference t p >t 
mean s.d. (= precision) 
Conductivity 39 5.0 0.04 1.79 0.09 
Gran plot alk. -6 39.2 0.06 1.24 0.23 
Double ept. alk. -7 38.7 0.15 1.69 0.11 
Suifate 260 18.7 0.01 0.95 0.36 
Color 13 3.3 0.01 0.27 0.79 
Hydrogen ion 12 11.6 0.06 0.66 0.51 





Double endpoint alkalinity vs. Gran plot alkalinity (unfrozen): 


(above) 0.07 0.90 0.38 





Sulfate was determined by the barium sulfate turbidimetric 
method (American Public Health Association et al. 1981) using a 
Hach model 2100A turbidimeter. Nine standards from 0-800 yeg/1 
SO4g were used to develop a calibration curve from which sample 
values could be determined. 


Chloride was determined using the mercuric. nitgate 
titrimetric method recommended by USEPA (Kopp and McKee 1979)°. 


Color was measured as an index of possible influence of 
organic acids, using the platinum-cobalt scale as represented by 
standard color discs in a Hellige Aqua Tester. 


> Haines and Akielaszek (1983) measured sulfate by the 
microthorin method (Fritz and Yamamura 1955) and chloride by 
the ferricyanide method (American Public Health Association et 
al. 1981). 
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On the advice of the auth-sr (James Kramer, McMaster 
University, personal communication), caicite saturation index 
(CSI) was calculatea from the modified equation given by Kramer 
(1981) rather than the general equation given by Kramer (1976). 


Mean amount and annual pH of precipitation at each site were 
determined from official state climatology maps and from the 
latest National Atmospheric Deposition Program (NADP) map 
summary, respectively. 


QUALITY ASSURANCE 


Analytical instruments were calibrated frequently against 
known standards and blanks. All analyses were performed in 
duplicate and statistical comparisons and evaluations were 
performed between the duplicates for each analysis. 


As a further check on analytical accuracy and data coding 
errors, we calculated total ionic balance and theoretical 
conductivity. All ion concentrations were converted _ to 
microequivalents per liter (yeq/l), anions and cations were 
summed, and the sum of anions was compared to the sum of cations 
for each sample. Stations with large ionic concentration 
imbalance (+ >100%) were removed from the data set before data 
were summarized. Theoretical conductivity was calculated by 
multiplying equivalent conductance of each ion from tables in 
Weast (1980) by the concentration of that ion for each sample and 
summing the results. The result was compared to measured 
conductivity for each sample as a check on the ionic balance; 
i.@., good conductivity agreement implies that all important 
ionic concentrations have been measured (Golterman et al. 1978). 


DATA ANALYSIS 


Results obtained from the above analyses were coded on a 
Standard form used by both the New England and Middle Atlantic 
Studies. Punched cards were produced and sent to Mr. Paul Rago 
at the U. S. Fish and Wildlife Service Great Lakes Fishery 
Laboratory in Ann Arbor, MI. Mr. Rago used the University of 
Michigan's MIDAS (Michigan Interactive Data Analysis System) to 
perform some of the more advanced statistical analyses and to 
produce the maps used in this’ report. Many of the basic 
correlations, regressions, and comparisons for this report were 
performed using a duplicate data set and the Statistical Analysis 











System (sas) ® as installed on the IBM computer complex at The 
Pennsylvania State University. Further error checking was done 
when the data were transformed for entry into the ACID 
(Acidification Chemistry Information Database) system at 
Brookhaven National Laboratory, as a contribution to the USEPA 
effort to determine sensitivity to acidification of waters 
throughout the United States (Omernik and Powers 1982). 


For lake and impoundment sites, results from the duplicate 
Surface samples were averaged together with results from the 
duplicate samples taken near bottom, to produce a single value 
for each measurement before further statistical analysis. for 
Streams, the duplicate mid-depth sample results were similarly 
averaged. 


In order to avoid undetected effects of unusual chemistry or 
possible contamination, all stations with color 30 units or 
greater, with double endpoint alkalinity greater than 1000 yeq/1l, 
with chloride greater than 500 ywegq/l, or with conductivity 
greater than 200 ywS/cm were eliminated from the data summaries 
which follow. These were in addition to those removed for ionic 
imbalance as described above. Many of the 60 stations so removed 
from the data set failed more than one of the above criteria. 
This process left 218 stations for further analysis. A list of 
the removed stations will be found in Appendix C. 


ee ee 





RESULTS AND DISCUSSION 


QUALITY ASSURANCE 


Comparison of the sums of anions and cations for each sample 
from the final data set yielded reasonable agreement (Figure 2). 
These calculations would not be expected to show perfect balance 
because not all the ions present were measured; anions such as 
Mitrate in particular may have been present in significant 
amounts in some waters. The regression line had an intercept of 
86.54 and a slope of 0.817 (R* = 0.621, p = 0.0001). # The 
intercept was significantly different from zero (p = 0.008) and, 
being positive, supports the hypothesis of “missing” anions. The 
Slope was not significantly different from 1.0 (p = 0.0001). The 
correlation was highly significant (r = 0.788, p = 0.0001). 


A further check was made by calculating the theoretical 
conductivity of each sample based on content of each major ion, 
and comparing these to the conductivities actually measured. The 
regression line (Figure 3) had an intercept of 3.260, slightly 
different from zero (p = 0.023), and a slope of 0.605, which was 
Significantly different from 1.0 (p = 0.0001) (R* = 0.847, p = 
0.0001). The correlation, however, was highly significant (r = 
0.920, p = 0.0001). Figure 3 also supports the hypothesis that 
not all ions present were measured, since in nearly every case 
the measured conductivity was greater than the theoretical based 
on those ions actually measured. 


We also performed an analysis of correlation and percent 
difference between the duplicate analyses for each sample. The 
results are presented in Table 4. 


These comparisons indicate that the data set is free of 
major analytical or transcription errors. 
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Figure 2. Sum of cations versus sum of anions for 211 sites, 
all samples averaged for each site. Least squares 
regression equation: sum of cations = 0.817 sum of 
anions + 86.54 (R’* = 0.621, p = 0.0001). 
Correlation coefficient (r) = 0.788, p = 0.001. 
(seven sites omitted due to missing values). 
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Figure 3: Theoretical versus measured conductivity for 211 
Sites, all samples averaged for each site. Least 
Squares regression equation: theoretical conductivity 
= 3.260 + 0.605 measured conductivity (R* = 0.847, 
p = 0.0001). Correlation coefficient (r) =0.920, 
p = 0.0001. (Seven sites omitted due to missing 
values). 
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Table 4. Comparison of Duplicate Chemical Analyses for 
all Samples 





Sample 1 vs Sample 2 





Ion or Number r mean standardized 
Measurement of pairs (p = 0.0001) % difference std. dev. 
( = precision) 





Hydrogen ion 325 0.974 1.4 57.6 
Alkalinity 318 0.967 3.2 65.1 
Conductivity 322 0.974 0.3 12.9 
Aluminum 333 0.845 17.8 55.0 
Calcium 333 0.999 ) OP 2.8 
Chloride 317 0.834 1.4 48.5 
Potassium 333 0.999 1.4 2.7 
Magnesium 333 0.999 1.2 2.6 
Manganese 333 0.999 42.9 81.2 
Sodium 333 0.999 1.6 19.2 
Sulfate 324 0.803 0.1 24.2 
Color 324 0.948 4.6 59.5 





SITE DISTRIBUTION 


We sampled 278 lakes and streams in the nine states. 
(Figure 1). The distribution of sites among states (Table 5) was 
roughly proportional to the areas of the states, but additional 
Stations were included where convenient (e.g., Delaware) or where 
waters were being sampled for other projects (e.g., 
Pennsylvania). Lakes and impoundments comprised 61% of the 
Sites, and streams comprised 39%. These percentages are somewhat 
different from those of the New England study (85% and 15% 
respectively) due to the scarcity of suitable lakes in most of 
Pennsylvania and some other areas. Stream order ranged from 1 to 
4 (N: 1=26, 2=44, 3=26, 4=14). Site elevations ranged from 2 m 
(MSL) to 1181 m, with a mean of 330 m. 


The distribution of sites by bedrock geology and soil 
classes present in each state is shown in Table 6. The 
distributions are not exactly proportional to each other or to 
the total area of each type in each state because suitable waters 
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Table 5. 


Number and Classification of Lakes and Streams Sampled in Each State 








Area of State 


No. (%) 


removed from 





State Total as ”, of total chemical 

(7) 9-state area analysis 

Impoundment summaries 
De laware 2 ») 21(7.6) 1.9 2(19) 
Maryland 17 0 24(8.4) 5.2 9(38) 
New Jersey 3 16 31(11.2) 3.9 8 (26) 
New York 8 28 46(16.5) 24.4 11(24) 
North Carolina 3 l 19(3.6) 6.5a 1(19) 
Pennsylvania 59 9 88 (31.7) 22.3 14(16) 
Tennessee 3 0 11(3.9) 5.2a 3(27) 
Virginia q 2 24(8.6) 19.64 4(17) 
West Virginia 6 l 23(8. 3) 11.9 8 (35) 
Total 119 278(190) 190.0 69 (22) 











“Areas for Tennessee and North Carolina were divided by four since only the Apnalachian parts of those 
states were sampled. 











Table 6. Distribution of Sampling Sites Among States by Bedrock and Soil Types 





Numbers and (%) of 
Stations on Soil Type 


Number and (%) of 


Number of Stations on Bedrock Type 
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_State _ Stations I II III VI NS SS1 SS2 S2 

Delaware 21 0 4 17 0 1 6 14 
(19) (81) (5) (28) (67) 

Maryland 24 0 15 5 4 8 19 5 1 

(62) (21) (17) (33) (42) (21) (4) 

New Jersey 31 5 24 0 2 12 9 10 N 
(16) (76) (8) (39) (29) (32) 

New York 46 7 32 2 5 31 11 2 2 

(15) (70) (4) (11) (68) (24) (4) (4) 

North Carolina 10 1 7 2 0 1 8 1 0 
(10) (70) (20) (10) (80) (10) 

Pennsylvania 88 0 51 25 12 46 21 15 6 

(58) (28) (14) (52) (24) (17) (7) 

Tennessee ll 0 3 ») 8 0 8 3 0 
(27) (73) (73) (27) 

Virginia 24 4 9 5 6 0 21 3 0 
(17) (37) (21) (25) (88) (12) 

West Virginia 23 0 10 12 l 6 16 1 0 
(44) (52) (4) (26) (79) (4) 

POrAL 278 17 155 68 38 105 110 54 19 

(6) (56) (24) (14) (38) (40) (19) (3) 





2) 











were seldom ideally located with respect to either factor. Sites 
were chosen primarily in proportion to bedrock and secondarily in 
proportion to soils. Most sites were on nonsensitive or slightly 
sensitive soils and/or on bedrock of low to medium acid 
neutralizing capacity. 


Distribution of sites by bedrock class is shown in Figure 4. 
As might be expected, sites in mountainous areas (generally in 
the northwestern half of the study area) are usually on bedrock 
classes I and II (most sensitive). On the other hand, there are 
frequent occurrences of sites on less-sensitive bedrock in 
mountainous areas. This emphasizes the point that generalized 
maps showing large sensitive (or nonsensitive) areas, e.g. Likens 
et al. (1979) and Comptroller General (1981), can be misleading. 


Distribution of sites by soil class is shown in Figure 5. 
No suitable waters on soil type Sl (“sensitive soils dominate the 
area") were available. 


CHEMICAL FACTORS 


Organic Acids and Color 





We measured color as an index of possible presence of 
Organic acids or other contaminants. Only eight of our 278 
waters had average color values of 30 or greater, and these were 
omitted from our data analysis. (Another 52 were omitted for 
other reasons.) For the remaining 218 stations, the mean color 
value was 6.04 units, and only six were above 20 units. 


Haines and Akielaszek (1983) derived the relationship TOC 
(mg/l) = 1.32 + 0.0613 color (Pt-Co units) from analysis of total 
Organic carbon (TOC) on a 3l-sample subset of their data. We 
used this equation for our samples and found no significant 
relationship between log)gTOC (R*’ = 0.00003, p = 0.980) or color 
(R? = 0.0001, p = 0.866) and pH. Similarly, there was no 
significant relationship between log TOC (R* = 0.0052, p = 0.276) 
or color (R* = 0.0028, p = 0.422) and alkalinity. 


We also examined the relationship between aluminum and 
Organic carbon, since it has been reported that aluminum forms 


strong organic complexes which could affect aluminum 
concentrations (Driscoll et al. 1980). The regression of 
aluminum on TOC for our samples was not significant (R’ = 0.0103, 
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Figure 4. Map showing distribution of sampling sites by 
bedrock classification (see Table 1 for description 


of classes). 
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Figure 5. Map showing distribution of sampling sites by soil 


Classification (see Table 1 for description 
of classes). 
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p = 0.126), nor was there a relationship between aluminum and 
color itself (R*’ = 0.0078, p = 0.184). 


pH and Alkalinity 





The distribution of waters by pH (surface values) is shown 
in Figure 6. Only in five of the nine states did we sample 
waters whose pH's were less than 5.5 at the time of sampling. 
Nearly half (13 of 31) of the waters with pH less than 5.5 were 
in the sandy, low-lying areas of New Jersey. Seven states had 
some waters with pH between 5.5 and 6.5. In all but one state, 
the largest group of waters was in the range pH 6.5-7.5. The 
exception was New Jersey, in which slightly more were below 5.5. 
Table 7 is a summary of the distribution of pH values by state. 


The distribution of alkalinity values (Figure 7) was similar 
to that of pH, although a somewhat larger portion fell in the 
more sensitive (lower alkalinity) class. Only Delaware had no 
waters with alkalinity below i100 yweq/l. There were more waters 
which appeared to be sensitive, or already becoming acidified, on 
the basis of alkalinity than on the basis of pH. This is not 
Surprising considering the temporal variation common in pH 
readings. Table 8 is a summary of alkalinity values by state. 


PH was correlated strongly (r >0.3) with alkalinity, 
conductivity and (negative) calcite saturation index, and weakly 
(r <0.3) with calcium, magnesium, and disturbance. Alkalinity 
was strongly correlated with calcium, pH, conductivity, and 
magnesium, and weakly with sulfate, chloride, disturbance, rain 
pH, and (negative) elevation. Neither pH nor alkalinity was 
correlated with aluminum, manganese, color, total organic carbon 
(as predicted from color), elevation, or size of waters. 


Haines and Akielaszek (1983) compared a number of regional 
surface water chemistry studies with respect to pH and 
alkalinity. Tables 9 and 10 are summaries of those studies 
compared with the present study; they suggest that waters with 
low pH and alkalinity are more likely to be found in areas 
receiving pronounced acidic precipitation. 


The overall pH of precipitation averaged 4.3 over our study 
area, but many locations within it averaged near 4.0. Our study 
waters appear generally to be less sensitive than most sites in 
other surveys of areas receiving acidic deposition, but are 
clearly more sensitive or, in some cases, more acidified, than 
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Figure 6. Map showing distribution of sampling sites by water 
DH. Data are for surface samples only and are means 
of duplicate analyses. 
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Table 7. Distribution of Average pH Values for Sampled 
Waters by State 





Percent in pH Range 








State n 

<§.5 5.5-6. 6.5-7.5 >7.5 

(n=20) (n=33) (n=112) (n=53) 
Delaware 19 0 ll 89 0 
Maryland 15 7 20 60 13 
New Jersey 23 40 30 30 0 
New York 35 3 l 51 35 
North Carolina 9 0 0 67 33 
Pennsylvania 74 ll 15 42 32 
Tennessee 8 0 0 75 25 
Virginia 20 0 10 50 40 
West Virginia 15 7 27 53 13 
All waters 218 9 15 51 25 








Table 8. Distribution of Average Alkalinity (peq/1) 
Values for Sampled Waters by State 





Percent in Alkalinity Range 














State n ee _ 
<100 100-200 200-400 >400 

(n=68) (n=36) (n=61) (n=53) 

Delaware 19 0 10 74 16 
Maryland 15 13 13 40 34 
New Jersey 23 57 17 9 17 
New York 35 26 20 23 31 
North Carolina 9 Ll ll 45 33 
Pennsylvania 74 43 16 1S 26 
Tennessee 8 318 0 25 37 
Virginia 20 l 20 50 20 
West Virginia 15 40 27 27 6 
All waters 218 31 17 28 24 
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Figure 7. Map showing distribution of sampling sites by 
alkalinity (yeq/l). Data are for surface 
samples only and are means of duplicate analyses. 
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DH of Surface Waters Reported by Regional Surveys 


(Modified from Haines and Akielaszek 1983). 





Location 


Number 


Percent 


in DH Range 





<5 


Reference 





Areas Where Precipitation Averages < pH 4.6 


Middle Atlantic 218 
New Encland 226 
West Sweden 314 
West Sweden 15 
Soutn Sweden 51 
South Norway 155 
South Norway 719 
Denmarx 14 
Scotland 72 
Nova Scotia 21 
Quebec 25 
Central Ontari 26 
LaClocke Mts. 152 
Sudbury, Ont. 150 


Adirondack Mts. 849 


28 


13 
25 


34 


iS 
30 


J 


2 
1 


' 


~ 


43 
27 
78 


72 
45 


This study 
Haines & 

Akielaszek 1983 
Almer et al. i974 
Dickson 1975 
Malmer 1975 
Wright et al. 1977 
Wright and 

Snekvik 1978 
Rebsdorf 1980 
Wright et al. 1980 
Watt et al. 1979 
Jones et al. 1980 
Scheider et al. 

1979 
Beamish and 

Harvey 1972 
Conroy et al. 1976 
Pfeiffer and 

Festa 1980 


Areas Where Precipitation Averages > pH 4.6 


North Norway 77 
N. W. Wisconsin 265 
North Minnesota 85 


0 


0 


13 


~ 


0 


87 


34 


100 


Wright and 
Gjessing 1976 
Lillie and Mason 

1980 
Glass and Loucks 
1980 














Table 10. Alkalinity of Surface Waters Reported by Regionai 
Surveys. (Alkalinity Values are in peq/1) 
(Modified from Haines and Akielaszek 1983). 





Percent in alkalinity range 
Location Number Reference 





<20 20-100 101-200 »>200 





Areas Where Precipitation Averages < pH 4.6 


Middle 218 15 17 17 51 This study 
Atlantic 
New England 226 23 18 12 47 Haines and 
Akielaszek 
(1983) 
South Norway 62 3 ll 2 84 Wright et al. 
1977 
Denmark 14 86 14 0 0 Rebsdorf 1980 
Nova Scotia 21 71 24 0 5 Watt et al. 
1979 
Central 26 12 73 15 0 Scheider et 
Ontario al. 1979 
Ontario 600 -- 164 32° 52 Zimmerman and 
Harvey 1979 
La Cloche 4 100 0 0 0 Beamish 1976 
Mountains 
Adirondack 692 41 25 18 16 Pfeiffer and 
Mountains Festa 1980 


Areas Where Precipitation Averages > pH 4.6 








Northwest 265 -- i5¢ 17 68 Lillie and 
Wisconsin Mason 1980 
North 85 0 22 26 52 Glass and 
Minnesota Loucks 1980 
a4 <50 yweq/l 


51-200 peq/l 
© <100 peq/l 











those receiving less acidic deposition. It is clear that many 
waters sensitive to acidification exist within the Middle 


Atlantic States. 


A map of acid deposition (Figure 8) shows relatively unifo 
loadings over the study area except for higher loadings near ¢t! 
common border of New York, New Jersey, and Pennsylvania. Thi 
border area, however, contained many waters with higher 
(Figure 6) and alkalinity (Figure 7) and low calcite Saturat ion 
index (Figure 9). These waters appear to be relatively resistant 
to acidification despite being located on relatively sensitiv 
bedrock types (Figure 2), perhaps because of nonsensitive and 
Slightly sensitive soils present in their watersheds (Appendix A) 
as noted by Gmur (Nicholas Gmur, Brookhaven National Laboratory, 
personal communication). 
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Calcite Saturation Index 








The calcite saturation index (CSI) has been proposed as an 
index of vulnerability to acidification by Conroy et al. (1974) 
and Kramer (1976). It is calculated as: 


CSI = p(Ca~*) + p(alkalinity) - p(H*) + DK 


- 


where: p(X) = 1990)9(X) 
p(K = 1.98 at 25 °C 
Ca = moles/l 


> 


alkalinity and H = equivalents,1 


Generally, higher values of CSI reflect lower ac} 
neutralizing capacity. We grouped our samples into categories of 
CSI less than 1, between 1 and 3, and greater than 3 (modified 
from Kramer 1976 and Glass and Loucks 1980). 


The distribution of stations by CSI is shown in Figure 9. 
It can be readily seen by comparison with Figures 7 and 8 that 
the areas where most waters had CSI less than 1 were in general 
those with the highest pH and/or alkalinity waters. Values of 
CSI iess than 1 represent waters that are saturated or nearly 
Saturated with respect to calcite and would not be susceptible to 
acidification from acidic precipitation. Values of 1-3 represent 
waters potentially susceptidle, and values greater than 3 
represent definite susceptibdil:.ty. Of the waters in our study 
11% were not susceptible, 57% were potentially | and 
32% were susceptibie, on the basis of CSI. CSI w 
>0.3) negatively with pH, alkalinity, calcium, c 
magnesium, as expected. There were weaker correlations with 











~ = 


10 


Figure 8. Map showing mean annual acid deposition from 
precipitation for period 1976-1979. (Units are 
mg/m’/yr hydrogen ion). From February 1981 
Interim Reports by work groups under U.S. - Canada 
Memorandum of Intent on Transboundary Air Pollution. 
Reprinted in Comptroller General (1981). 
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Figure 9. Map showing distribution of sampling sites by 
calcite saturation index. Data are for surface 
samples only and are means of duplicate analyses. 
Higher values of CSI indicate lower acid 
neutralizing capacity. 
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elevation (positive), disturbance, sulfate, and amount of annual 
precipitation (all negative). Table 11 is a summary of CSI 
distribution by state. 





Table 11. Percent Distribution of Calcite Saturation 
Index by State 





Calcite Saturation Index 











State n 
<1 1-3 >3 

Delaware 19 0 89 ll 
Maryland 15 13 67 20 
New Jersey 23 4 65 31 
New York 35 23 49 28 
North Carolina 9 12 44 44 
Pennsylvania 74 12 54 34 
Tennessee 8 0 75 25 
Virginia 20 5 60 35 
West Virginia 15 7 27 66 
Whole Study 218 ll 57 32 





Haines and Akielaszek (1983) determined that 59% of the 
waters they surveyed in New England could be classified as 
susceptible on the basis of CSI. They reviewed several other 
Studies (Kramer 1981, Glass and Loucks 1980, and Zimmerman and 
Harvey 1979) and concluded that CSI does not provide a predictive 
ability any better than that of pH and alkalinity measurements 
alone. In our study, 49% of the waters sampled were susceptible 
on the basis of alkalinity less than 200 yeq/l, and 32% were 
Susceptible on the basis of CSI greater than 3. Furthermore, CSI 
showed a strong negative correlation with alkalinity (r = -0.657, 
p = 0.0001) (Figure 10). Therefore, we see no advantage in the 
use of CSI as a predictor of vulnerability of Middle Atlantic 
States waters over the use of alkalinity, except that it does 
relate alkalinity and vulnerability to calcium concentration. 


32 


¥\ 











Log Alkalinity lyeq-I ‘| 


Figur’ 








Calcite Saturation Index 


Relation of log of alkalinity to calcite saturation 
index for 218 sampled waters (one point hidden). 
Data are means of duplicate analyses and means of 
surface and bottom samples where taken. Least 
squares regression equation: log alkalinity = 

2.367 - 0.154 CSI (R* = 0.038, p = 0.004). 
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Conductivity and Ionic Composition 





Conductivity (= specific conductance) is an easily measured, 
widely used parameter which is directly related to the total 
ionic content of a water sample. It thus has potential as a 
Simple index of sensitivity to acidification. However, in some 
cases, factors such as proximity to the seacoast (a source of 
airborne salts) and differences in ion activities may limit its 
usefulness as a predictive device. 


The 218 waters included in this analysis had a range of 
conductivity from 17 to 199 uS/cm, and averaged 70. Conductivity 
was positively correlated (r >0.03) with alkalinity, calcium, 
chloride, magnesium, sodium, and sulfate, and negatively 
correlated with elevation, pH, CSI, Secchi disc transparency, and 
disturbance. The relationship between conductivity and 
alkalinity (Figure 1l) was statistically significant (R* = 0.502, 
p = 0.0001) but a few waters with relatively high conductivity 
were low in alkalinity. This may be due to conductivity of 
elevated hydrogen ion (Weast 1980) or unusual mineralogy. There 
waS no correlation between conductivity and distance from the 
seacoast. Table 12 is a comparison of the relative concentration 
of major ions in our waters with those of the New England study 
(Haines and Akielaszek 1983) and with world averages. 





Table 12. Percent Relative Composition of Major Ions in Sampled 
Waters Compared with New England (Haines and Akielaszek 
1983) and World Averages (three world averages from 
other authors combined from Henriksen 1980) 





New Englaid This study 
Ions World average 








DH<S pH5-6 pH>é PH<5S pH5-6 pH>6 








Cations 
Ca+Mg 88.2 79.5 79.9 91.9 62.8 75.6 80.9 
Na+K 11.8 20.5 20.1 9.0 37.2 24.4 19.1 
Anions 
HCO3 82.7 0 15.4 74.8 0 12.5 56.9 
SO4 17.3 100 84.6 25.2 100 87.5 43.1 
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Figure ll. Relation of alkalinity to conductivity for 218 
Sampled waters (2 points hidden). Data are means 
of duplicate analyses and of all depths sampled. 
Least squares regression equation: conductivity = 
0.115 alkalinity + 39.225 (R* = 0.502, p = 0. 0001). 








As noted in the New England investigation, ionic composition 
of the surface waters varied markedly with pH (Figure 12). 
However, all pH classes of our waters had higher mean ionic 
content than the corresponding New England class, and the 
differences between total means were not as great. Total ionic 
content was highest at high pH and declined at successively lower 
pH. 


Other similarities between our data and that from New 
England are evident: bicarbonate disappeared below pH 5, as 
expected, and sulfate remained at approximately the same 
concentration while making up a greater proportion of total 
anions as pH decreased. The cations calcium and magnesium 
declined markedly with pH, but sodium and potassium showed no 
such trend. This is probably related to the sources of these 
ions; calcium and magnesium are assumed to originate primarily 
from the watershed and are related to acid neutralizing capacity, 
but sodium and potassium are assumed to originate primarily from 
other sources and are not involved with buffering reactions. 
Aluminum was approximately equal in the pH >6 and 5-6 classes, 
but more than doubled below pH 5. This is consistent with the 
fact that aluminum is more soluble at low pH and thus is 
frequently mobilized from watersheds under acidic conditions. 


Chloride 





Chloride was measured aS a potential indicator of 
disturbance by human activity and of influence of sea salts. 
Haines and Akielaszek (1983) found that chloride concentrations 
were a poor indicator of human disturbance for New England 
waters. They also found a strong inverse relationship between 
chloride and distance to the seacoast, and hence used chloride 
concentration as a correction for the ionic contribution of 
Marine aerosols. This correction technique was suggested by 
Wright and Gjessing (1976), who estimated seawater sulfate by 
assuming that all the chloride comes from seawater and that the 
amount of sulfate attributable to sea salts is proportional to 
the sulfate/chloride ratio in seawater. The seawater 
contribution of the other ions can be estimated in similar 
fashion. Wright and Gjessing (1976) suggested that the 
correction was important for areas near the coast but relatively 
minor for inland lakes. Also, this technique involves’ the 
assumption that no sulfate is retained or released by soils. 


We believe that the seaspray contribution to lakes in 
eastern North America is less than the contribution to lakes in 
Scandinavia (where Wright and Gjessing worked) due to our 


36 











700 


600 


500 


400 


300 


200 


100 


Figure 12. 
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<5 5-5 >6 


Mean ionic composition of sampled waters classified 
by pH range. For each water, data are means of 
duplicate analyses and of all depths sampled. 

(N = number of waters in each range). 
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prevailing westerly winds. Our data do not show a strong 
correlation between chloride or conductivity and distance to the 
Seacoast. We used the ratio of 55% Cl to 7.7% SOqg, giving a 14% 
correction to be subtracted for sulfate from marine aerosols; 
i.@., nonmarine SOqg = totai SOy - 0.14 Cl. Aiso, nonmarine Ca = 
total Ca - 0.021 Cl and nonmarine Mg = total Mg - 0.067 Cl 
(Sverdrup et al. 1942). In the following sections we have 
selectively included, and so labeled, both corrected 
("nonmarine") and uncorrected data. 


Chloride concentrations in our final set of waters ranged 
from 1 to 473 weq/l, with a mean of 124. Chloride was positively 
(p = 0.0001) correlated with disturbance, alkalinity, calcium, 
magnesium, conductivity, sodium, and sulfate, and negatively 
correlated with elevation. The latter is to be expected if 
higher elevation, headwater streams tend to be both less 
disturbed and more sensitive to acidification. Of the ten waters 
from our final set that were highest in chloride, none were near 
the seacoast and all but two had high disturbance codes. One 
case with low disturbance was the largest lake sampled (684 ha), 
which is a limited-access water supply but lies near an area of 
commercial salt production (with possible underwater exposure to 
salt beds); the other was a stream. 


Sulfate 





Sulfate values ranged from i8 to 932 yueq/l (averages of 
duplicate analyses at all depths). Mean sulfate concentration 
for the entire study was 227 yweq/l. Sulfate was positively (p = 
0.001) correlated wih alkalinity, conductivity, calcium, 
chloride, magnesium, manganese, and sodium; but was less stronaly 
correlated with disturbance (p = 0.07). It had a weak inverse 
correlation with calcite saturation index. All these relations 
are as expected, and there is no evidence of any significant 
contribution from marine aerosols. When corrected for possible 
marine aerosol contribution in the manner of the New England 
study (Haines and Akielaszek 1983), our sulfate values were 
independent of chloride. 


The concentrations of sulfate found in the Middle Atlantic 
States are higher than those reported for some other areas where 
precipitation is acidic, and much higher than those for areas 
with nonacidic precipitation (Table 13). As in the New England 
Study, neither total nor nonmarine sulfate was significantly 
correlated with hydrogen ion concentration (r = 0.026, 0.029 
respectively, p = Q.5). Like the New England results, ou 
calcium and calcium+¢magnesium values were weakiy correlated wit 
Sulfate (r = 0.384, 0.464 respectively, p = 0.0001). 
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Table 13. Concentrations of Sulfate (yeq/l) Reported 
from Surface Waters by Regional Surveys 
Total Sulfate Nonmarine Sulfate 
Region n mean s.d. mean s.d. Reference 





Areas where precipitation averages < pH 4.6 


Middle 

Atlantic 218 227 161 211 158 
New 

England 226 138 51 127 48 
South 

Norway 26 99 33 92 - 
West 

Sweden 6 200 70 180 - 
Scotland 72 178 9} 148 76 
Nova Scotia 16 152 28 117 25 
Central 

Ontario 15 198 97 - - 
Sudbury, 

Ontario 4 804 292 800 290 
La Cloche 

Mountains 4 287 42 285 41 


Areas where precipitation averages > pH 


This study 
Haines and 
Akielaszek 
1983 
Gjessing et 
al. 1976 
Grahn et al. 
1974 
Wright and 
Gjessing 1976 
Watt et al. 
1979 
Dillon et al. 
1978 
Scheider et 
al. 1975 


Beamish 1976 


4.6 





Central Gjessing et 
Norway $2 65 34 63 - al. 1976 
Western Armstrong and 
Ontario 7 63 = 58 - Schindler 

1971 
39 
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Although all stations in the coastal areas of our study (New 
Jersey, Delaware, and eastern Maryland) showed sulfate levels 
above 100 yweq/l (Figure 13), this was also the case in several 
noncoastal areas (western New York, western Pennsylvania, western 
West Virginia; in general the Allegheny Plateau'. A man of acid 
deposition over the region (Figure 8) shows tnat except for a 
small area of high Geposition in the New York-New 
Jersey-Pennsylvania boundary area, deposition over the study area 
is relatively uniform. The variations in sulfute are probably 
more dependent on geolcgical factors such as gypsum or pyrite 
deposits than on acidic deposition. 


Calcium and Magnesium 





We discuss calcium and magnesium tocether due to their 
Similar origin and behavior. The sum of calcium and magnesium 
was correlated with a.kalinity (r = 9.640, p = 0.0001), Dut not 
as strongly as in the New England Study (r = 0.95) (Haines and 
Akielaszek 1983) nor as in Henriksen (1980) (r = 0.98). This may 
be due to the influence of some higher alkalinity waters on our 
data. The regression equation of our data (Figure 14) is 
alkalinity = 0.42 (Ca*Mg) + 75.9 (both in ueq/l) (R*% = 0.409, p = 
0.0001). Thi intercept (+75. st. is consideczably higher than 
those reported in the above-mentioned studies, reflecting the 
generally hicgher level of alkalinity in the Middle Atlantic area. 


When the streams and impoundments are removed from our data, 
leaving only “natural” lakes, the intercept drops to 23, still 
not as low as those of the New England study (-55) (Haines and 
Akielaszek 1983) and Henriksen's (1980) study of data from areas 
relatively unaffected by acid precipitation. The slope of our 
Gata is lower than for either of the above studies, suggesting a 
lesser role of calcium+magnesium in the higher alkalinity wate:is 
of our area. This is similar i:o the observations of Almer et al. 
(1978), that the alkalinity:(Ca+Mg) relationship was near 1:1 
where sulfur deposition was low, but less than 1:1 where 
deposition was high. Calcium concentrations for Middle Atlantic 
States waters are compared with those found in studies of other 
areas in Table 14. 


Aluminum and Manganese 





Aluminum was not significantly currelated with pH or 
hydrogen ion concentration (r = -0.166, p >0.01 and r = 0.143, p 
>0.03, respectively). As expected, high aluminum concentrations 
were generally found in low pH waters, but there were enough 
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Figure 13. 
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Relation of alkalinity to calcium + magnesium for 
218 sampled waters (one point hidden). Regression 
equation is alkalinity = 0.418(Ca+Mg) + 75.99 

(R*? = 0.409, p = 0.0001). 
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Table 14. Concentrations of Calcium Reported from Surface 
Waters by Regional Surveys (Modified from Haines 
and Akielaszek (1983)) 





Calcium content, yeq/l 
Region n Reference 





Mean S.D. Range 





Areas where precipitation averages < pH 4.6) 


Middle Atlantic 218 279 242 21-1279 This study 
- lakes only 45 321 255 43-1086 
- streams only 85 269 248 25-1279 
- impoundments 88 266 229 21-1031 
only 
New England 226 270 395 16-2300 Haines and 
Akielaszek 
1983 
South Norway 26 57 33 - Gjessing et 
al. 1976 
Nova Scotia 16 124 98 42-453 Watt et al. 
1979 
Scotland 72 239 371 14-1815 Wright and 
Henriksen 
1980 


Areas where precipitation averages > pH 4.6) 


World average - 239 - - Henriksen 
1980 
Central Norway 52 139 110 - Gjessing et 
al. 1976 
Western Ontario 40 80 - 30-225 Armstrong and 
Schindler 
1971 





43 


3EST COPY AVAILABLE 














exceptions to prevent a showing of statistical significance. The 
regression equation for our waters is log Al (peq/l) = 2.83-0.23 
pH (R* = 0.048, p = 0.0011). Both this intercept and slope are 
lower than those reported by other regional surveys (see Table 12 
of Haines and Akieiaszek 1983). In general, however, our 
aluminum values were higher than those found in New England, 
ranging from 0 to 395 yweq/l, with a mean of 46.8. When the 
"natural" lakes or the streams are separated from the rest of our 
data, the relationship for each alone is only slightly stronger 
(R? = 0.123 and 0.114 respectively, p = 0.001). Waters with high 
aluminum concentrations were widely distributed throughout our 
Study area (Figure 15). Only six of 218 stations (2.8%) had 
aluminum concentrations over 200 peq/l, and only one of these was 
at a pH below 5.5.° This is the pH level below which elevated 
aluminum concentrations are considered toxic to some species of 
Fish (Schofield and Trojnar 1980). 


Although Haines and  Akielaszek (1983) felt that the 
Similarity of relations between pH and aluminum concentration in 
Studies reviewed by them would enable the prediction of aluminum 
concentration from pH in any region, this would apparently not be 
the case in the Middle Atlantic States. This may be due to the 
generally higher ionic content, leading to more interferences, 
and/or to the generally greater complexity of bedrock and soil 
types in our area. 


Manganese and aluminum are chemically similar, and our 
results for the manganese-pH relation were similar to those for 


aluminum. The regression equation is logjgMn = 0.11-0.045 pH, 
(R? = 0.0017, p = 90.03), a weaker relationship than that for 
aluminum. Haines and Akielaszek (1983) found a similar lack of 
correlati_.. bet ‘een manganese and pH in New England waters. 


Manganese vaiues in our study ranged from 0 to 52 yveq/l, with a 
mean of only 3.2, much lower than the aluminum concentrations. 


Neither aluminum nor manganese was strongly correlated with 


alkalinity (r = -0.166, 0.181 and p = 90.0140, 0.0076 
respectively), mor were aluminum and manganese significantly 
correlated with each other (r = -0.055, p = 0.416) despite their 


chemical similarity. 


~~ —S— eee Se ee ee ee or or 


7 Note that in their discussion in the corresponding section, 
Haines and Akielaszek (1983) used different units of 
concentration. 
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Map showing distribution of sampling sites by 
aluminum concentration in weq/l. Data are 
means of duplicate analyses and of surface 


and bottom samples where taken. 
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Precipitation Chemistry 





As shown in Figure 8, our study area had some variation in 
acid loading, but differences in bedrock, soils, and other 
influences masked any relationship which might exist between acid 
loading and water chemistry. A loading index was derived from 
the product of mean annual weighted precipitation pH (National 
Atmospheric Deposition Program 1984) and annual rainfall (from 
State climatological summaries for the year of sampling) (values 
listed in Appendix A). Both this index and its logarithm were 
tested against pH, alkalinity, elevation, conductivity, aluminum, 
calcium, chloride, color, and sulfate. In all cases the 
correlation was insignificant (r <0.02). 


PHYSICAL FACTORS 


Elevation 





Neither pH nor alkalinity was significantly correlated with 
elevation, although 24 of the 32 stations at elevations greater 
than 590 m had alkalinities below 200 yeq/l, the normal “cutoff 
point" for sensitivity. Neither was there any significant 
relationship between elevation and bedrock sensitivity, soil 
sensitivity or disturbance. It was noted, however, that all 1l 
Stations on bedrock class I (very low acid neutralizing capacity) 
occurred in the lower half of the elevation range (< 590 M). 
This was somewhat contrury to our expectations, and appeared to 
be due to their location on the lower slopes of the poorly 
buffered rock areas. Most such areas in the Middle Atlantic 
States occur on hill and mountain tops where permanent water 
bodies are rare. There were strong negative correlations between 
elevation and chloride and conductivity (r = <-0.512, -0.451 
respectively, p = 0.0001). 


As did Haines and Akielaszek (1983) in New England, we found 
acidic waters at low elevations, and waters with pH less than 5.5 
(some below 5.0) at elevations from near sea level to over 1,000 
m. We also found 21 waters with pH greater than 6.5 at 
elevations above 590 m; such waters seldom occurred in the New 
England study. Although it is certainly true that higher 
elevation areas uSually have thin soil and sensitive bedrock, 
other factors (perhaps disturbance) apparently prevent these from 
being dominant influences on the chemistry of Middle Atlantic 
waters. This is in contrast to the situation in Norway (Wright 
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and Snekvik 1978) where acidic lakes are rare at lower 
elevations, apparently due to the presence of easily weathered 
sedimentary bedrock. 


Size 





Neither pH nor alkalinity was significantly correlated with 
lake or pond area or stream width. Most of the natural lakes 
sampled were less than 50 ha in area (mean 47) and had pH below 
7. Most of the impoundments were also less than 50 ha in area 
(mean 36), and their mean pH's tended to be lower than those of 
the “natural” lakes. Most small lakes had low alkalinity, but 
there were numerous small, higher alkalinity lakes and a few 
large, low alkalinity lakes. 


Stream width is related oniy imprecisely to discharge or 
other measures which might be expected to reflect’ acid 
neutralizing capacity. Only six of our 85 stream stations had 
widths greater than 10 m; five of these six had pH greater than 
7.5, and their alkalinities ranged from 227 to 842 yeq/l. It is 
generally true that wider streams have larger drainage areas, so 
that greater buffering by natural processes in the watershed and, 


possibly, human influence are expected. A high negative 
correlation has been found between stream width and acidification 
in Pennsylvania streams (Fred Johnson, Pennsylvania Fish 


Commission, Harrisburg; unpublished). 


Hydrology Type 





Both natural lakes and impoundments were separated into two 
types - seepage lakes (with no inlet) and drainage lakes (with 
inlet). In general, seepage lakes occurred at higher elevations, 
were smaller, and were less disturbed, but most chemical factors 
were not significantly different between the two types (Table 
15). 


The chemical factors measured, except alkalinity and 
calcium, were higher in the drainage lakes, probably reflecting 
leaching from the upstream watershed as well as the immediate 
lake basin. The differences, however, were not statistically 
Significant. Drainage lakes were generally lower in elevation, 
larger, and had higher disturbance values, as might be expected. 
These differences are probably responsible for the higher ionic 
concentrations observed in the drainage lakes. Eilers et al. 
(1983) found that in northern Wisconsin, lake hydrology was a 
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Table 15. Chemical and Physical Factors for Lakes Classified by 
Hydrology Type. "“p" is the probability that means 
are drawn from the same population (t-test of one- 
way analysis of variance). (No significant 
differences found) 





Hydrology type 

















Factor Drainage Seepage 
(n=119) (n=14) 
Mean $.5d. Mean S.D. p 
pH, units (from H+ conc.) 5.6 0.9 6.7 0.6 0.61 
alkalinity, peq/l 256 227 360 303 0.70 
conductivity, p»pS/cm 71 40 65 38 0.22 
calcite sat. index, units 2.7 1.3 2.5 1.2 0.65 
calcium, yedq/l 281 242 317 215 0.49 
magnesium, yeq/l 182 131 152 124 0.28 
sodium, ypeq/l 67 53 48 64 0.49 
potassium, yeq/l 47 27 31 27 0.25 
aluminum, yeq/l 48 56 14 12 0.04 
manganese, yeq/l 4 8 5 5 0.69 
chloride, ypeq/l 141 131 50 55 0.27 
sulfate, yeq/l 233 173 208 185 0.21 
color, Pt-Co units 7.0 5.6 5.3 2.8 0.50 
Secchi disc, m 1.8 Lea 1.9 1.3 0.56 
disturbance, units 3.0 0.9 2.6 0.9 0.06 
elevation, m MSL 285 238 428 209 0.12 
area, ha 42 84 24 31 0.22 





Significant factor in water chemistry; but Haines and Akielaszek 
(1983) did not find any significant effect of hydrology in New 
England. Differences between the chemistry of lakes and streams 
(Table 16) were significant only for calcite saturation index, 
manganese, color, and disturbance. 


Haines and Akielaszek (1983) also found that streams were 
higher in pH than were lakes, but the only other ionic 
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Table 16. Chemical and Physical Factors for Lakes and Streams. 
"p" is the probability that means are drawn from 
the same population (t-test of one-way analysis of 

















variance) 
Lakes Streams 
(n=133) (n=85) 
Factor p 
Mean S.D Mean S.D. 

pH, units (based on H* conc.) 5.6 0.9 5.7 1.0 0.06 
alkalinity, yeq/l 267 237 263 240 0.73 
conductivity, yuS/cm 70 39 69 38 0.42 
calcite sat. index, 2.6 Les 203 1.3 0.09 
calcium, yeq/l 284 238 269 248 0.58 
magnesium, yeq/l 179 148 166 130 0.72 
sodium, peq/l 66 54 68 66 0.83 
potassium, yeq/l 46 27 47 28 0.64 
aluminum, yeq/1l 44 54 51 61 0.53 
manganese, yeq/1l 4 8 l 3 0.00 
chloride, yeq/l 134 128 112 110 0.90 
Sulfate, peq/l 230 174 222 140 0.44 
color, Pt-Co units 6.8 5.4 4.2 3.8 0.00 
Secchi disc, m 1.8 1.2 NA NA NA 
disturbance, units 3.0 0.9 2.6 0.9 0.01 
elevation, m MSL 300 238 343 211 0.33 
area, ha; or width, 40 80 5.7 5.1 NA 
concentration differences between streams and lakes their 


Study were in sodium and potassium rather than in manganese. 


It 


seems unlikely that any of these differences are of importance in 
thus we have pooled 
data from all waters for many of our analyses. 


controlling vulnerability to acidification, 


Sharpe et al. 


found 


that 


pH, 


alkalinity, 
If true 


dissolved aluminum were functions of stream discharge. 


for our streams, 


this would 


indicate 


that 


and 


the conditions we 


observed would be influenced by the flow existing at the time we 


sampled. 
lesser extent. 


Unfortunately, 


neither 


resources to do periodic sampling were available to us. 
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This problem probably exists for lakes as well but to a 
discharge data nor 


the 











Stream Order 





Average width increased with stream order, but no 
Significant (at the 0.05 level) differences were found amonc 
streams of different orders for any of the characteristics 
measured, except disturbance. There was, however, a general 
increasing trend in content of calcium and chloride, and a 
decrease in calcite saturation index and ievation with 
increasing stream order (Table 17). This is the expected 
pattern, due to the larger, more varied watersheds of higher 
order streams. Johnson (unpublished) examined 1,812 coldwater 
stream sections in Pennsylvania and found a high negative 
correlation between stream order and three measures of 
acidification (pH < 6, alkalinity < 5 ppm, alkalinity < 10 ppm). 





Table 17. Chemical and Physical Factors for Streams by 
Stream Order. Means with same following letter 
are significantly different (one-way analysis of 
variance with Scheffe's comparison, p <0.05) 





Stream order 











Factor l 2 3 4 
(n=18) (n=35) (n=21) (n=11) 
pH, units 5.5 5.6 6.1 5.9 
alkalinity, ypeq/l 279 215 280 356 
conductivity, yuS/cm 66 61 69 95 
calcite sat. index, units 2.2 2.5 2.3 1.5 
calcium, yeq/l 213 261 261 407 
magnesium, peq/l 144 175 153 197 
sodium, yveq/l 65 82 53 54 
potassium, ypeq/l 49 59 43 34 
aluminum, yeq/1l 53 64 38 26 
manganese, yeq/l l l 2 l 
chloride, ypeq/1l 70 104 107 207 
sulfate, peq/l 236 223 213 213 
color, Pt-Co units 2.6 4.6 4.4 5.2 
disturbance, units 2.4 2.3a 2.9 3. °a 
elevation, m MSL 405 374 335 164 
width, m 2.9 4.3 6.1 14.1 
50 








Bedrock Classes 





Analysis of chemical results compared among bedrock 
sensitivity classes showed significant differences only in 
alkalinity and potassium (Table 18). This is similar to the 
results from the New England study except that they found 
Significant differences in pH, calcium, conductivity, and calcite 
Saturation index as weli. There was no steady change in any 
factor from bedrock class I (very low acid neutralizing capacity) 
to IV (very high acid neutralizing capacity) as might have been 
expected. This may be related to the fact noted above (under 
Elevation) that the class I bedrock stations were not at the 
highest elevations (all were less than 550 Mm). Due to the 
relatively small number of sampling sites on class I bedrock, it 
is not advisable to draw conclusions from these data concerning 
the influence of interactions between bedrock and elevation on 
water chemistry. The pH of waters located over bedrock of class 
II or III varied over the whole range 4-8 (Figure 16). Only 
waters with pH above 6.25 were found over bedrock of class IV, 
and, unexpectedly, over bedrock of class I. 


Examination of the relationship between bedrock and 
alkalinity levels (Figure 17) yields a similar impression. 
Although class I bedrock has the lowest acid neutralizing 
Capacity, no waters with very low alkalinity (< 75 wpeq/1) 
occurred in class I bedrock areas. Most waters in class II 
bedrock areas had low alkalinities, and waters in classes III and 
IV bedrock areas had alkalinities over nearly the full range 
considered. Thus even though many waters did reflect the 
relative acid neutralizing capacity of their underlying bedrock, 
the relationship was not consistent. Haines and Akielaszek 
(1983) obtained similar results in New England. They pointed out 
that the bedrock classes proposed by Hendrey et al. (1980) are 
related to chemical factors associated with acid neutralizing 
capacity and acidity but not to other factors which might 
influence water chemistry. 


Because of the lack of significant differences between 
waters on bedrock classes I and II and similarly between III and 
IV, classes I and II were pooled and classes III and IV were 
pooled for further statistical analysis. Class I-II was termed 
“sensitive” and class III-IV was termed "“nonsensitive."” Also, 
waters with alkalinity of 200 yeq/l or less were termed 
"sensitive" and those higher in alkalinity were termed 
“nonsensitive."” A similar grouping was done for pH, with the 
division placed at pH 6.75. We tested the effect of bedrock on 


51 


re 


OV UP AWE 











—— 





Frequency 







HI 








10 


Lnormt lh 


45 50 55 60 — 70 7.5 80 
p 





Figure 16. Frequency distribution of waters by bedrock 
geology in pH intervals. pH values listed 
are midpoints of ranges of 0.5 unit. 
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Figure 17. Frequency distribution of waters by bedrock 
geology in alkalinity intervals. Alkalinity 
values listed are midpoints of ranges of 
50 weg/l. 
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Tabie 18. Chemical and Physical Factors for Waters Classified 
by Bedrock Geology. Bedrock classes are: I = very 
low acic neutralizing capacity; II = moderately low 
acid neutralizing capacity; III = moderately high 
acid neutralizing capacity; IV = very high acid 
neutralizing capacity (Hendrey et al. 1980). Means 
with the same fcllowing letter are significantly 
different (one-way analysis of variance with Scheffe's 
comparison, p = 0.05) (Three waters in class V omitted) 





Bedrock class 








Factor I II III IV 
(n-11) (n=127) (n=54) (n=23) 

pH, units (from H* conc.) 6.9 5.5 6.1 7.2 
alkalinity, yeq/l 424 229a 263 415a 
conductivity, uS/cm 75 65 78 78 
calcite sat. index, units 2.2 2.5 2.6 2.3 
calcium, yweq/l 304 122 231 370 
magnesium, yweq/!l 206 175 165 165 
sodium, yeq/l 109 70 62 43 
potassium, yweq/l 69a 49 45 34a 
aluminum, yeq/1 81 49 43 30 
manganese, yweq/l 6 3 3 2 
chloride, yeq/l 127 122 154 75 
sulfate, yweq/l 210 228 233 238 
color, Pt~-Co units 5.4 5.6 7.3 4.2 
disturbance, units 3.5 2.7 3.0 3.0 
eievation, m 260 342 241 346 
Secchi disc (lakes only), m 1.4 1.7 1.§ 2.7 





pH and alkalinity using several statistical measures based on 
contingency tables. The most useful of these were the phi 
measure (Sokal and Rohlf 1973) and the odds ratio (Bliss 1967; 
Fleiss 1973). The contingency table format allows us to ask the 
following questions: 
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1. Are the numbers of observations within categories any 
Gifferent than what would be expected due to chance 
alone? 


2. Is there any significant association between the bedrock 
Or soil class and the pK and/or alkalinity of the 
overlying waters? 


3. What are the odds of observing a sensitive water, given 
that the bedrock or soil type is sensitive? 


The general format of such tables is: 


Water category 











Predictor Sensitive Non-sensitive Sum 
Sensitive Nl Nn}2 nN}. 
Non-sensitive N 2} nz? n>. 
Sum Ne} n.? n.. 


If the variables are fully independent, then the observed cell 
frequencies will be the products of their respective row and 
column marginal probabilities. Chi-square (X*) is the usual 
Statistic for testing the null hypothesis of independence in such 
tables. The association measure (phi) is also used to evaluate 
predictive usefulness; in general associations less than 0.35 are 
not considered particularly useful as predictors (Sokal and Rohlf 
1973). Another, more appropriate measure is known as the odds 
satio (Fleiss 1973). If A represents a sensitive bedrock or soil 
Class and B represents a sensitive lake type, then the odds ratio 
(w) gives the odds, or chance, of observing B when A is present 
relative to the chance of observing B when A is absent. 


These analyses for ““irock are summarized in Tables 19 and 
20. The results agree » ti the earlier indications (Table 17) 
that bedrock is only a marginally u*eful predictor of the 
sensitivity of a water to acidificstion (based on pH and 
alkalinity). The odds of observing a sensitive lake in a 
sensitive bedrock area are only 1.1 to 1 and 1.9 to 1 for pH and 
alkalinity, respectively. Oniy in the case of alkalinity does 
the chi-square value indicate that the ratio is significantly 
different from i. 
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Table 19. Relations Between Bedrock Sensitivity and Water 
Sensitivity Based on pH 














pH 
Bedrock Sum 
Classes 
<6.75 >6.75 
xX? = 0.057 
{(p = 0.812) 
I é& Ii! 54 84 138 
III & IV 30 30 80 phi = 0.016 
w = 1.071 
Sum 84 134 218 s.e.(w) = 0. 








Table 20. Relations Between Bedrock Sensitivity and Water 
Sensitivity Based on Alkalinity 














Alkalinity 
Bedrock Sum 
Classes 
<200 >200 
x‘ = §.277 
(p = 0.0216) 
I & I! 74 64 138 
III & IV 30 50 80 phi = 0.156 
w = 1.927 
Sum 104 114 218 s.e.(w) = 0.553 





Soil Classes 





Results of chemical analyses of sampled waters grouped by 
SOil class (i.@., SOii Cation exchange capacity of McFee (1980) 
as defined in Table 1) are presented in Table 21. Most factors 
exhibited a significant difference among soil classes, whereas in 
the New England study only calcium, chloride, sulfate, sodium, 
and alkalinity were significantly different. Although we have 
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placed class SSl second in rank in Table 21, waters there exhibit 
the lowest mean calcium concentration and alkalinity; they also 
have the hichest mean calcite saturation index and lowest mean 
conductivity (along with those on class S2) and the lowest mean 
pH. On the other hand, our means for physical factors suggest 
that the order of sensitivity might be S2>SS1>SS2>NS, as we have 
presented it. 





Table 21. Chemical and Physical Factors for Waters Classified 
by Soil Sensitivity. Soil Classes are: 
Sl= sensitive; S2= sensitive <50%; SSl=slightly 
sensitive; SS2=slightly sensitive <50%, NS= non- 
sensitive. Means with the same following letter 
are significantly different (one-way analysis of 
variance with Scheffe's comparison, p = 0.05) 
(No waters were sampled on class Sl soil) 





Soil class 








Factor 
S2 SSsl SS2 NS 
(n=39) (n=86) (n=46) (n=47) 
pH, units (from H*conc.) 6.2 5.5 5.6 5.7 
alkalinity, yeq/l 242 228 324 297 
conductivity, uS/cm 57a 57ab 96a 78b 
calcite sat. index, units 2.9 2.7 2.2 2.2 
calcium, yeq/1l 240 193ab 268b 379a 
magnesium, yeq/l 106ab 142cd 24lad 223bc 
sodium, ypeq/l 34ac 58b 90ab 86c 
potassium, yeq/l 34ab 50a 5lb 48 
aluminum, ypeq/1l 3la 63ab 44 32b 
manganese, yed/l 3 3 5 3 
chloride, yeq/l 78ab 83cd 199bd l66ac 
sulfate, peq/l 183a 194b 325abc 229c 
color, Pt-Co units 4.8 5.1 6.9 7.0 
disturbance, units 2.5a 2.8 3.1la 2.9 
elevation, m MSL 434ab 33l1lc 217bc 293a 
Secchi disc, m (lakes) 2.7a 1.8 1.5a hel 











McFee's (1980) class definitions (Table 1) present problems 
of interpretation when we attempt to assign individual waters to 
a class on the map. For example, a lake may clearly be located 
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in an area labeled "S2" on the McFee map, Dut the full definition 
of “S2" is: “sensitive soils are significant, but cover less than 
50% of the area.“ Is our lake thus on sensitive soil or not? 
Only a local, detailed survey of the soils around a water body 
can provide such information; even then the picture may be 
clouded by the presence of several adjacent soil types. Unknown 
local soil effects may mask relationships between soil class and 
water chemistry. Also, the maps at the scale available to us 
were not sufficiently precise to associate a given water with an 
exact soil class with certainty in every case. To a lesser 
extent, the same problem may exist for the bedrock classes. It 
Should be noted that the original intent of McFee's (1980) study 
was not that the classes be used as we have, but to as.ess the 
proportions of total area in the eastern U.S. likely to be 
sensitive on the basis of soil cation exchange capacity. 


The frequency distribution of pH of waters by soil classes 
(Figure 18) is quite similar to that for bedrock. The greatest 
numbers were in the interval pH 6.76-7.25, and classes SSl, SS2, 
and NS occurred in nearly all pH intervals. Class $2 occurred 
only above pH 4.75. Most, but not ali, waters on soil class NS 
had pH of 6.25 or greater. This distribution was similar to that 
founc by Haines and Akielaszek (1983) except that they had no 
waters on soil class S2. The relationship beteen soil classes 
and alkalinity (Figure 19) was not as skewed as it was for New 
England, except that the largest number of waters were in the 
highest and lowest alkalinity intervals (>375 and <25 yeq/l). 
All four soil types were represented in nearly every alkalinity 
interval. 


Contingency tables were also developed for soil and water 
sensitivity on the bases of pH and alkalinity. As before, pH was 
divided above and below pH 6.75, and alkalinity above and below 
200 wpeq/l. For this analysis, soil classes S2, SSl, and SS2 were 
combined and termed "sensitive," and class NS alone was termed 
"nonsensitive.” Results are given for pH and alkalinity in 
Tables 22 and 23, respectively. The contingency tables for soil 
class are not significantly different from what would be expected 
by chance alone. The odds ratios (1.2:1 for pH and 1.1:1 for 
alkalinity) both indicate that lake sensitivity is not strongly 
related to soil type. 


Haines and Akielaszek (1983) found that the odds ratios for 
both bedrock and soil class provided good predictions of 
sensitivity, although bedrock was much better than soil. Kaplan 
et al. (1981) concluded that bedrock influenced water chemistry 
not directly but rather through soils. However, they used 
different classifications of bedrock and soil from ours, and 
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Figure 18. Frequency distribution of pH of waters by soil 


classes. pH figures given are midpoints of 
intervals of 0.5 pH. 
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Figure 19. Frequency distribution of alkalinity of waters 
by soil classes. Alkalinities are midpoints of 


intervals of 50 yeq/1l. 
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Table 22. Relations Between Soil Sensitivity and Water 
Sensitivity Based on pH 














pH 
Soil 
Classes 
<6.75 >6.75 Sum 
x? = 0.409 
(p = 0.5225) 
S2 + SSl1 + SS2 20 27 171 
NS 64 107 47 phi = 0.043 
w = 1.238 
Sum 84 134 218 s.e.(w) = 0.414 








Table 23. Relations Between Soil Sensitivity and Water 
Sensitivity Based on Alkalinity 





Alkalinity, yueq/l 














Soil 
Classes 
<200 >200 Sum 
_ KX? = 0.036 
(p = 0.8489) 
SZ + $Sl + SS2 23 24 47 
NS 81 90 171 phi = 0.013 
ow = 1.065 
Sum 104 114 218 s.e.(w) = 0.35 





studied primarily lowland sites where the soils would be expected 
to be thicker than those in our study and the companion New 
England work. Thin soils such as are found at higher elevations 
would be expected to allow a greater influence of bedrock on 
water chemistry. Sensitive soils and sensitive bedrock may occur 
together and sometimes act synergistically in influencing (or 
failing to influence) the acidity of precipitation runoff, the 
chemistry of receiving waters, or both. Thus we again conclude 
that, at the presently available level of detail, knowledge of 
only the bedrock and soil conditions is not a sufficient basis 
for predicting the likelihood of acidification of a given water 
body, and that the published maps showing sensitive areas based 
On bedrock alone can be misleading. 


61 





Human Disturbance 





The degree of huméen disturbance in the watershed was 
Significantly related to. elevation, pH, conductivity, and 
chloride at the 95% confidence level (Table 24). This was 
Similar to the New Englanc study results. 

The least disturbed waters were at the highest elevations, 
and the most disturbed were at the lowest. For lakes, the Secchi 
disc transparency was higher in the least disturbed lakes and 
lower in those most disturbed, as expected. Most of the chemical 
factors closely related to acidification - alkalinity, calcite 
Saturation index, calcium, magnesium, and sulfate - showed 
consistent trends, and pH was higher in classes 3 and 4 than in 
classes 1 and 2. Aluminum, potassium, sodium, color, and 
manganese were not relatec to disturbance in any consistent way. 
Unlike the New England results, color in our study was not a 
potential index of disturbance. Chloride, however, may have some 
use in this regard. 

Human activities such as construction, waste disposal, 
farming, logging, and road travel result in exposure of soils and 
rock to increased leaching of chemical constituents, and also may 
increase the amount of surface runoff. This is, no doubt, the 
cause of these observed relations. Unlike Haines and Akielaszek 
(1983), we found no statistical evidence of a relationship 
between bedrock class and disturbance (Table 18), but soil 
Classes S2 and SS2 were significantly related to ‘@isturbance 
(Table 21). The overall correlation between soil class and 
disturbance, however, was low (r = 0.134, p = 0.048). 


MULTIVARIATE ANALYSES 


In any study with as many variables as the present one, it 
is desirable to explore muitivariate analysis techniques in order 
to clarify the importance of particular variables. We performed 
several such techniques of which the results of two are 
presented. 


Principal Components Analysis 





Principal components analysis is a multivariate technique 
for examining relationships among several quantitative variables. 
We used it to identify those chemical factors which seem to have 
the most influence on the overall water chemistry, and thus might 
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Table 24. Chemical and Phys.cal Factors for Waters by 
Degree o: duman “isturbance (for definitions of 
classes 1-4, see Table 2). Means with the same 
following letter are significantly different (One- 
way analysis of variance with Scheffe's comparison, 











p = 0.905) 
Disturbance class 
Factor 
l 2 3 4 
{n=21) (nm=51) (n=88) (n=58) 
pH, units (from H* conc.) 5.2a 5.6 5.9a 5.8 
alkalinity, yweq/l 144 222 273 334 
conductivity, wS/cm 46a 61b 71 84ab 
calcite sat. index, units 2.9 2.6 2.6 2.1 
calcium, yeq/l 261 248 291 311 
magnesium, yveq/l 122 151 172 215 
sodium, peq/l 56 52 70 77 
potassium, yeq/l 45 46 47 49 
aluminum, ypeg/1l 40 52 38 58 
manganese, ypeq/l 3 2 4 3 
chloride, yeq/l 57a 78b 131 174ab 
sulfate, peq/l 174 209 240 244 
color, Pt-Cc units 4.9 5.4 6.2 6.1 
Secchi disc (lakes only), m 1.8 2e3 1.8 1.5 
elevation, m MSL 535a 399b 315 214ab 





have value as predictors of vulnerability to acidification. 
Results of this analysis for the major chemical measurements of 
our study are presented in Table 25. 


The first principal component accounts for the largest 
portion of the variance and is dominated by magnesium, calcium, 
chloride, sodium, and alkalinity. Sodium and chloride are not 
involved with acid neutralizing capacity or acidity, but may 
reflect disturbance or the presence of salt deposits in some 
watersheds. Alkalinity, calcium, and magnesium are closely 
related to acid neutralizing capacity. The second component is 
dominated by alkalinity and calcium, with aluminum and potassium 
showing an inverse influence. The third component is dominated 
by aluminum and potassium, and inversely by hydrogen ion, the 
Significance of which is unclear. Results for the New England 
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Table 25. Principal Components Analysis of Chemical Factors 





Principal Components 











Variable 

First Second Third 
Cumulative percent 
of varianc® explained 25 40 48 
by componerit 
Hydrogen ion -0.07 -0.31 -0.40 
Alkalinity 0.29 0.30 0.15 
Aluminum -0.01 -0.33 0.38 
Calcium 0.35 0.17 -0.20 
Chloride 0.32 -0.05 -0.09 
Potassium 0.12 -0.43 0.35 
Magnesium 0.39 0.03 -0.08 
Sodium 0.30 -0.32 0.12 
Sulfate 0.26 -0.02 -0.15 





data were similar except that sodium and chloride were important 
only in the second component, and hydrogen ion (as pH) was 
important in the first component. 


From these results, it appears again that alkalinity is 
potentially useful as a predictor, while calcium, magnesium, and 
sulfate are also strongly influential on the chemistry of the 
sampled waters. The principal components technique is relatively 
weak for this data since the first three components account for 
only 48 per cent of the variance. 


Canonical Correlation Analysis 





Canonical correlation is a technique for examining the 
relationship between two sets of variables by testing the 
hypothesis that each canonical correlation and ail those smaller 
are actually equal to zero in the population. Given two sets of 
variatles, a linear combination is found from each set, such that 
the correlation between the two combinations, called canonical 
variables, is maximized. This is called the first canonical 
correlation. The coefficients of the combinations are called 
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canonical weights, and can be used to determine which variables 
have the most influence on the variance of the data set. Next a 
second set of canonical variables, uncorrelated with the first 
pair, is found to produce the second highest correlation 
coefficient. This process continues until the number of canonical 
variables equals the number of variables in the smaller group, 
but usually not all the pairs are used (SAS Institute 1982). 


We compared chemical measurements to physical factors from 
our data using canonical correlation analysis. Only the first 
three canonical correlations showed reasonably high probability 
of being different from zero, and are presented in Table 26. The 
results show good correlations (r, = 0.70, rg = 0.35, rz = 0.31) 
between physical and chemical factors, all significant at the 
0.05 level. Soil class and disturbance carry the highest weights 
among the physical factors, while chloride, sodium, and 
alkalinity dominate the chemical factors. This is the first case 
in our data analyses wherein soil appeared to be a stronger 
influence than bedrock. 


We performed a canonical redundancy analysis (van den 
Wollenberg 1977) which examines how weli the original variables 
can be predicted from the canonical variables. This analysis 
shows that neither of the first pair of canonical variables is a 
good overall predictor of the opposite set of variables, the 
proportions of variance explained being only 0.103 (chemical) and 
0.135 (physical). The first three canonical variables together 
explain only 0.134 and 0.175 of the chemical and physical factor 
variances, respectively. The squared multiple correlations 
calculated as part of the redundancy analysis indicate that the 
first canonical variable of the physical factors has some 
predictive ability for chloride (0.313) and sodium (0.271) but 
little or none for the other ions. The first canonical variable 
of the chemical factors tas some predictive power for elevation 
(0.392), but little else. 


We conclude that canonical correlation analysis is not a 
Strong predictive tool for use with our data, but in general does 
not contradict the conclusions drawn from our other methods. 


HISTORICAL COMPARISONS 


At least one set of historical data was found for 160 of the 
278 waters in our sampling program. For the following analysis, 
we used unscreened data from our original set of 278 waters, but 
only a few comparisons involved waters not in the _ reduced 
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Table 26. Canonical Correlation Analysis of Standardized 


Chemical and 


— » | = 
hysical Factors 





Canonical correlation 








Factor 
First Second Third 
Canonical correlations 0.70 0.35 0.31 
(p (corr = 0)) 
Physical factor weights: 
size 0.09 0.40 0.10 
elevation -0.77 0.2 0.24 
bedrock class -0.19 0.74 0.14 
$0Oil class 0.27 0.17 0.94 
Gisturbance 0.22 0.55 -0.50 
Chemical measurement weights: 
hydrogen ion 0.22 -0.46 0.13 
alkalinity 0.30 0.49 “0.73 
aluminum 0.17 0.09 -0.45 
calcium -0.25 0.11 0.60 
chloride 0.65 0.23 -0.36 
potassium -0.20 -0.52 -0.06 
magnesium 0.07 -0.29 -0.64 
manganese 0.17 -0.40 -0.2 
sodium 0.64 -0.00 -0.09 
Sulfate 0.19 0.33 -0.04 
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(screened) set. In several cases multiple historical data points 
exist, bringing the total to about 350 sets of historical values. 
Most o” these contained pH data, but many were lacking alkalinity 
data, and a few had alkalinity but not pH. The oldest data found 
were from 1930. A listing of all historical data obtained for 
each station is given in Appendix 8. 


pH Comparisons 


We compared our pH data with the oldest available pH data 
for each station (Figure 20). The low coefficient of 
determination (R* = 0.288) indicates that the changes in pH over 
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Recent pH 





Figure 20. 








Oldest Historical pH 


Recent (this study) pH value versus oldest value 
available for 36 locations. Diagonal line 
indicates equal pH. Waters with pH > 8, alkalinity 
> 250 weq/l, or time span not greater than 2 

years have been excluded. Numbers represent time 


between observations grouped as follows: 1 = 2-5 
years; 2 = 6-10 years; 3 = 11-15 years; 4 = 16-20 
years; 5 = >20 years. 
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time are not consistent among sites. Of the 160 waters, 62 (393%) 
decreased in pH and 98 (61%) remained the same or increased. 
These are nearly the opposite of the comparable figures or New 
England waters as reported by Haines and Akielaszek (1983). As 
pointed out earlier by Arnold et al. (1960), there are a number 
of factors, mostly related to man's activity, that can be 
expected to cause increases in the pH of waters over time, but 
few that cause a decrease in pH other than inputs of acid from 
precipitation or obvious pollution sources such as new coal mine. 
in the watershed. In general, the greater incidence of human 
activity in the Middle Atlantic States compared with New England 
is probably an important reason for the higher number of cases in 
which pH or alkalinity increased. 


The length of time between th Oldest “historical” pH 
measurement and the measurement done by us at the same location 
varied from 1-50 years, with a mean of 17 years. We plotted 
change in pH against elapsed time in years for each site and 
found no significant trend (R*’ = 0.0003, p = 0.73). #£=This 
analysis included the entire set of waters for which historical 
Gata were available, and indicates that there is no consistency 
in rate of change, if any, among the sampled waters. Many of 
these lakes have high acid neutralizing capacity and would not be 
expected to show much change in pH, even with significant inputs 
of acidity. 


As any experienced field worker knows, measurement of 
accurate pH in the field, especially on poorly-buffered waters, 
is a task having limited and inconsistent success regardless of 
the method used. The historical data available to us sometimes 
specified the general method (electrometric or colorimetric) used 
and sometimes did not. Such information is not a sufficient 
basis for judging either the accuracy or the precision of the 
data. Furthermore, changing limnological conditions over a few 
hours or days, or between seasons, can cause significant changes 
in pH and thus render historical pH comparisons meaningless 
unless they involve numerous data points on the same water or 
large, consistent changes. Our experience leads us to conclude 
that the above analyses of changes and trends in the available 
data cannot be used to draw conclusions about overall trends 
because the accuracy of the historical data is unknown. 


Alkalinity Comparisons 








A more instructive approach is to lock at changes in acid 
Neutralizing capacity as represented by alkalinity. For these 
analyses, fixed endpoint alkalinity (by titration to pH 4.5) was 
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used for Ou “recent” measurement, Since considerable 
investigation DY ourselves and other workers has indicated that 
almost all alkalinity data before approximately 1975 were 
produced by fixed endpoint methods. Kramer and Tessier (1982) 
and Haines and Akielaszek (1983) presented reviews of the 
literature and considerations involved in historical alkalinity 
comparisons. The fixed endpoint method often overestimates true 
alkalinity (Kramer and Tessier 1982). we compared the fixed and 


double endpoint alkalinities of our samples and found that they 
were highly correlated (r = 0.9984, p <0.01), with the regression 
equation being double endpoint = 0.377 fixed endpoint - 25.3 (R*= 
0.997, p = 0.0001). Thus’ the ixed endpoint historical 


alkalinity could be co d for the overestimation by 
Subtracting 25.3 yweq/l. This correction was virtually identical 
when the data were separated as natural lakes, streams, or 
impoundments. Haines and Akielaszex (1983), using a factor of 32 
peq/l, compared corrected historical data to their double 
endpoint data for New England waters; whereas we measured the 
fixed endpoint aikalinity and compared it directly with the 
historical data. Kramer and Tessier (1982) suggested that the 
true correction factor to be subtracted from alkalinity data 
obtained by methyl orange titration is 81 peq/l. We were unable 
to determine exactly which of the historical data we obtained 
were done by methyl orange and which by electrometric titration, 
nor the exact endpoints used in mos” cases. Thus the results of 
these comparisons should be considered approximate. 


The mean historical aikalinity for all 95 available 
comparisons was 697 yeq/l, anc our mean “recent" alkalinity for 
the same waters was 489 yeq/l, or a decrease of 208 yveq/l (20%). 
As was done for pH, we plotted the change in alkalinity against 
the time in years between observations for each site and found no 
Significant relationship (R* = 1.0001, p = 0.88). This indicates 
that the rate of decrease in alkalinity varies considerably among 
the sampled waters. Only 26 ‘23%) waters showed an increase or 
no change in alkalinity, while 87 (77%) decreased. Haines and 
Akielaszek (1983) foui1d corresponding values of 30% and 70%, 
respectively, in New England. They reviewed several other 
Studies which all showed decreases in mean alkalinity of waters 
in Sweden and the United States over the same time period as our 
data. Johnson (1984) examined data for 41 Pennsylvania streams 
Surveyed by Pennsylvania Fish Commission personnel using the same 
techniques over a time span of about 15 years and showed that 92% 
decreased in alkalinity, from an average of 400 yeq/l to an 
average of 180 yeqg/1. 


In a number of articles on the “acid rain problem" (e. g., 
see papers by Hodanbosi, Frohliger, Smith and others in Kostik 
1980; also Curtis 1980 and Hunton and Williams 1983) the 
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Figure 21. Recent (this study) fixed endpoint alkalinity value 


versus Oldest available for 19 locations. Diagonal 
Line indicates equal alkalinity (no change). Waters 
with time span not greater than 2 years have been 
excluded. Numbers represent time between 
observations grouped as follows: 1 2-5 years, 

2 = 6-10 years, 3 = 11-15 years, 4 16-20 years, 

5 = > 20 years. 
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assertion has been made that the problem is not real because few 
waters have shown a major change in pH. However, because of the 
bicarbonate buffering system, little charse in pH is expected 
over a wide range of hydrogen ion lcading. Thus alkalinity 
changes are a better index of acidification thar are pH changes. 
Also, alkalinity is a better index of sensitivity to 
acidification at a given site, or to compare sensitivity between 
sites. Schofield (1982) determined -he change of pH between 
historic and recent data for lakes in the Adirondack Mountains of 
New York, and compared it with recent calcium concentrations. He 
calculated that at a level of 100-125 yweq/l calcium, the 
bicarbonate buffering capacity would be lost and pH would 
decrease. This agreed well with the changes which he observed. 
We did a similar comparison and observed some similar results, 
but we also observed some increases in pH below 125 ueq/l 
calcium, and a few decreases at higher calcium levels. These 
exceptions were all cases where the time span of the historical 
comparison was quite long (10-20 or more years). The exceptions 
may also support a hypothesis that has been suggested by Coombs 
(Marjorie Coombs, Florida Department of Environmental Regulation, 
personal communication) and others - that as acid neutralizing 
Capacity declines to critical levels, pH and other water 
chemistry measurements become highly variable among successive 
observations at a site. Only long-term periodic monitoring of 
Such sites can test this hypothesis. 


Although alkalinity lS, in general, a more Stable 
characteristic of a natural water than is pH, it can vary 
somewhat in time ard space wichin a water body, especially a 
poorly buffered one. Thus, as in the case of pH, we feel that 
only large, consistent changes based on several data points for 
the same water are useful to suggest acidification based on 
historical data. Even data showing these characteristics can 
only be evaluated if the exact methodclogy and titration 
endpoints used are known so that appropriate correction factors 
can be applied as discussed above. These cautions and the small 
number of alkalinity comparisons available to us suggest that the 
data shown in Figure 21 should be considered only a suggestion 
that a decrease in alkalinity may have occurred in many waters of 
the Middle Atlantic States. 


APPLICATION OF PROPOSED MODELS OF ACIDIFICATION 


Several models and ciassification systems based on wacer 


chemistry have been proposed to determine and/cr predict the 
degree of acidification of waters or their potential 
susceptibility to acidification. We applied several of these tc 


our data as discussed below. 











Sensitivity Classification 





Altshuller and McBean (1979) proposed a classification by 
groupings under four common water chemistry measurements: 
alkalinity, bicarbonate concentration, calcium corcentration, and 
conductivity. Since in our study bicarbonate was not determined 
separately from alkalinity we have not considered it here. We 
have also converted mg/l to yuweq/l for consistency (Table 27). 





Table 27. Sensitivity Classes of Altshuller and McBean (1979) 
(Percent of Waters in the Present Study Occurring 
in Each Shown in Parentheses) 














Alkalinity Calcium Conductivity All three 
Sensitivity factors 
(% of total) 
(weq/l) (uS/cm) 
High 0-200 0-200 0-30 
(48%) (513%) (9%) (93%) 
Moderate 200-400 200-400 31-70 
(28%) (22%) (47%) (3%) 
Low >400 >400 >70 
(24%) (26%) (44%) (16%) 





Although the agreement between distributions of sensitive 
waters based on alkalinity anc on calcium is fairly good, there 
is poor agreement between these distributions and that based on 
conductivity. A larger proportion of waters fall into the higher 
conductivity classes, probably due to the presence in many waters 
of elements such as _ chloride, sodium, and potassium that 
contribute to conductivity but not to acid neutralizing capacity. 
When the data are sorted in this way, only 9% of all waters meet 
all three criteria for high sensitivity, and only 3% meet all 
three criteria for moderate sensitivity. These percentages are 
much smaller than those of sensitivity determined by other 
Criteria discussed earlier. The classes based on alkalinity and 
calcium, however, agree well with our earlier ccnclusions 
regarcing sensitivity on the basis of alkalinity and calcite 
Saturation index. As shown earlier, the calc.te saturation index 
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is highly correlated with alkalinity and also may serve as a 
predictor of sensitivity, but it is more complex and has no clear 
advantage other than to link pH and calcium with alkaJinit 


Calcium-pH Model 





Henriksen (1979) proposed a model involving the relationship 
between pH and calcium content. He derived an empirical curve 
for this relation: waters falling above the curve are classified 
as acidified (i.e., they have a lower pH than would be expected 
on the basis of their calcium content, due to depletion of their 
acid neutralizing capacity). When this curve is applied to our 
data (Figure 22) for waters with nonmarine calcium of 300 yeq/1 
Or less, we find that 21 of 75 waters (28%) fall above the curve. 


As pointed out Dy Haines and Akielaszek (1983), the 
relationship used by Henriksen presumes that nonmarine calcium 
content of surface waters is not affected by acidic deposition 


Their linear regression of alkalinity against the sum of 
nonmarine calcium and maqnesium was alkalinity = -53 +¢ 
0.93(Ca+Mqg) (R* = 0.948), which agreed well with Henriksen's 
value of alkalinity = <-14 + 0.93(Ca+Mg), except for a lower 


intercept. The corresponding regression equation for our data 
(Figure 14) is alkalinity = 76 + 0.42(Ca+Mg) (R* = 0.409, p = 
0.0001), a higher intercept but a lower slope. The positive 
intercept for our data indicates that the calcium and magnesiu 
ions were generally balanced by bicarbonate, and its magnitude is 
probably due to the larger proportion of other anions such as 
chloride in our waters. 


In the New England study anions other than sulfate and 
bicarbonate appeared to be negligible as shown by ionic balance 
data which did not include chloride. Our ionic balance (Figure 
12) indicates a substantial chloride component, but chloride 
probably does not enter importantly into buffering reactions. 


Haines and Akielaszek (1983) presented a review of fr 
studies concerning the effect of acidic deposition on calci 
magnesium, or total cation, content of waters. c 
Suggests that sulfate increases are halanced by equivale 
bicarbonate decreases so iong as bicarbonate iS availabie, and 
that when bicarbonate is absent, sulfate increases are bDalancec 
by increases in cations such as aluminum, magnesium, and calcium. 
The lack of historical data on oncentrations of these ions 


+ 


Cc 
prevents us from testing this hypothesis directly. Figure 12 
shows that of these ions, only aluminum seemed to increase in 
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Figure 22. Calcium-pH model (Henriksen plot) applied to 
Middle Atlantic States lakes having calcium 
concentrations less than 300 yveq/l. Waters 
falling above the curve are considered acidified 
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concentration at low pH values where bicarbonate was absent. 
However, the data of Ficure 12 are means, and trends in 
individual waters over time are obscured, especially considering 
the generally higher total ionic content of the higher pH waters. 


Sulfate-pH Model 





Henriksen (1979, 1980) also proposed that acidification 
could be related to the nonmarine sulfate content of surface 
waters. He constructed a nomograph relating pH to calcium plus 
magnesium and sulfate content of waters. We applied the 
regressions of this nomograph to our lake daca and determined its 
success in predicting lake pH (Table 28, Figure 23). The 
predictive ability of the model was generally poor on data from 
the Middle Atlantic States, except for those lakes in the range 
of pH 4.7-5.3. The small number of our lakes in the “acidic" 
range makes it difficult to evaluate that part of the model. 
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Table 28. Prediction of pH of 82 Middle Atlantic Lakes Having 
Sums Of Nonmarine Calcium Plus Magnesium Less than 
500 peq/l by Henriksen's Sulfate-pH Model 
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pH range 
<4.7 4.7-5.3 >§.3 
Number actualiy in pH range 6 3 73 
Number predicted to be in pH range 8 22 50 
Number predicted correctly * l 3 48 
Percent predicted correctly * 17% 100% 66% 





* i.e., portion of those actually in the pH range that fall 
within the predicted isopleths on the plot 





Haines and Akielaszek (1983) performed a similar analysis 
and also reported poor predictive abiiity although the 
distribution of their results was somewhat different. They 


suggested that the poor predictive ability might be due to 
correlation of sulfate with pH in their data, or to an extraneous 
source of sulfate. They improved the model's predictive success 
somewhat by applying an arbitrary reduction factor to tne sulfate 
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Figure 23. Sulfate-pH model (Henriksen nomograph) applied to 
82 Middle Atlantic States lakes having calcium 
plus magnesium concentrations less than 500 yeq/l. 
Diagonal lines indicate Henriksen's (1980) 


regressions for waters of pH 4.6-4.8 ("4.7") and 
5.2-5.4 ("5.3"). Symbols indicate pH of waters 
sampled in this study: #® = pH < 4.7, () = pH 4.7- 


5.3, @= pH > 5.3. 








values. We did not experiment with the model in this way because 
our sulfate data were not correlated with pH and because, as 
Stated earlier, we found no evidence for a marine source of 
Sulfate in our data. The use of the “nonmarine” correction for 
Sulfate has already reduced our values below the actual 
measurements. 


Cation Denudation Model 





Thompson (1982) proposed a model of acidification based on 
the idea of “cation denudation” or the export of cations from a 
watershed in runoff. it is similar to the sulfate-pH model 
(Henriksen nomograph) except that the sum of cations is used 
rather than only calcium plus magnesium. With this model pH can 
be predicted from the carbonate buffering system of natural 
waters. 





Table 29. Prediction of pH of 74 Middle Atlantic Waters Having 
Sums of Nonmarine Cations Less Than 600 yeq/1 and 
Nonmarine Sulfate Less Than 500 yeq/l by Thompson's 
Cation Denudation Model 








DH range 

















at, 7 . 
lumber actually in 


pH range 7 10 3 Ll 37 
Number predicted to be 

in pH range 2 14 30 15 ll 
Number predicted 

correctly * 0 2 5 4 8 
Percent precicted 

-orrectiy * 0 20 56 36 57 

* ..@., portion of those actually in the pk range that 
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Figure 24. 
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Thompson cation denudation model applied to 74 
Middle Atlantic States waters with sums of 
cations and sulfate concentrations less than 
600 peg/i. Numbers iudicate observed pH 
ranges (sampie values averaged over depths and 
replicates): 1 4.0-5.0, Z 5.4°76.2, 3 
6.3-6.5, 4 6.6-6.8, 5 > 6.8 
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samples. The predictive success was poor, ranging from 0-57%. 
In the New England study the fit to this model was good at lower 
pH values but deteriorated as pH of samples increased. This was 
mot true of our data; in fact the agreement was better at higher 
PH values. There was no other obvious pattern to the agreement, 
and little improvement would be derived from reducing the sulfate 
values. 


As pointed out by Haines and Akielaszek (1983) this model 
depends on the relationship between excess sulfate and 
acidification, as does the Henriksen pH-sulfate model. It can 
thus be used to predict pH changes for given sulfate loads if the 
relationship between sulfate in surface waters and sulfate in 
atmospheric deposition is known. Unfortunately this is not yet 
the case. 


The cation denudation model can also be related to 
acidification which has already occurred. All points falling 
below the pH 5.1 line have other nonmarine cations in excess of 
bicarbonate; they may thus be classified as acidified. Of the 74 
waters in our survey which are included in Figure 24, only 2 (3%) 
fall in the “acidified” category on this basis. This is 
considerably --ss than the 28% found by the calcium-pH model, and 
the 39% which decreased in pH based on historical data; or, of 
course, the 77% which decreased in alkalinity. Thus it appears 
that this mode! is not particularly useful for studying 
acidification in waters of the Middle Atlantic states. 








6. 


CONCLUSIONS 


About 45% of the unpolluted, relatively undeveloped waters 
sampled in the nine Middle Atlantic States are vulnerable 
to acidification on the basis of alkalinity being less 
than 200 yweq/l, or calcium content less than 200 yueq/1l. 
Differences in water chemistry between lakes and streams, 
as groups, were small. 


Alkalinity is probably the most useful and relatively 
accurate predictor of sensitivity to acidification. 
Caicite saturation index is also a valid predictor, but 
has no advantage over alkalinity. 


Comparisons with historical data on the same waters indicated 
no significant overall change in pH but some decrease in 
the alkalinity of many waters. The accuracy of the 
historical data cannot be evaluated, however, so no firm 
conclusions regarding change can be drawn from it. 


The calcium-pH pilot of Henriksen is a fairly good predictor 
of lake pH values as related to calcium content. It 
predicted that about 28% of the sampled lakes» were 
acidiried. 


Neither bedrock nor soil class is a particularly strong 
predictor of sensitivity to acidification. Bedrock class 
is somewhat related to alkalinity of Middle Atlantic 
States waters, but soil type is not, at least when 
determined from rock and soil maps presently available. 
Thus some published maps showing geographic areas which 
are sensitive to acid precipitation, but which are based 
on bedrock alone, may be misleading. 


There was only one case of a potentially toxic aluminum/pH 
combination among the sampled waters. 
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RECOMMENDATIONS 


The principal unanswered questions raised in this report 
could best be addressed by the collection of water chemistry data 
on a regularly repeated basis, at least once per season, from a 
selected group of the waters sampled for this study. Waters with 
low levels of calcium, alkalinity, and total ionic content should 
be emphasized. At least a majority of the waters should have 
data on their past water chemistry available, in addition to that 
collected for this study. Also, concurrent data on the status of 
fish populations and their food chains in the same waters should 
be collected in order to document the effect, or lack thereof, of 
observed water chemistry trends on the aquatic ecosystem. 
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APPENDIX A. DATA OF THIS STUDY 


Data in the head block for each station are arranged as 
follows by line (see text for details): 


1. station number - water body name - county - state 

2. latitude (degrees, minutes, seconds) - longitude - USGS 
7.5-minute quadrangle name - elevation (m) 

3. date of sample (year, month, day) - Secchi disc 
transparency (m) (lakes only, streams = ***) 
- approximate area (ha) of pond/lake, or width (m) of 
stream - bedrock class (see Table 1) - soil class 
(see Table 1) 

4. disturbance code (see Table 2) - annual mean pH of 
precipitation at site - mean annual amount of 
precipitation at site (mm) 


The labeled columns are as follows: 


deep 


D depth of sample (m) 

T temperature at sample depth (°C) 

FLD PH pH measured in the field 

LAB PH pH measured in the laboratory 

FEP ALK fixed endpoint alkalinity as 
calcium carbonate (ypeq/1) 

DEP ALK double endpoint alkalinity as 
calcium carbonate (ypeq/]) 

COND conductivity in ywS/cm 

Al aluminum (total) in yveq/l 

Ca calcium in ypeq/l 

Cl chloride in peq/l 


tne second 


CLR color in Pt-Co units 

K potassium in peq/l 

Mg magnesium in yveq/1l 

Mn manganese in yeq/1 

Na sodium in peqg/l 

SO4 sulfate in peq/1l 

most cases where there are samples listed from 

decths, the second depth is the near-bottom sample taken near 
deepest part of the or pond. (In a few cases our samplinsc 
equipment would not bottom; these were generally over 


Sample was taken at the maximum 


possible but not necessarily near the bottom). 
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olf 7/22 1.J8 T6ha BR2 SS 
DiST+ RAINpH: 4.3 111888 


0 23 6.5 6.45 129 69 29 8 110 31 18 19 55 2.4 7 106 
2 23 666 6.99 1848 83 35 9 106 73.—=C:«dN#9 20 54 2.5 7 223 
2195 ALDER LAKE ULSTER NY 
cr 420259 744051 AR@NA 6742 
a B14 T1722 «1.58 16ha BB2 SS 
— DIST2 RAINpH: 4.3 116808 
omy 0 23 7.6 9.57 253 213 40 0 253 7 6 9 65 1.2 4 89 
= 4 17 6.64 6.47 287 245 67 0 242 26 6 19 66 4.9 & 111 
—_ © 2196 SUNHILL POND DELAWARE wY 
= “i 422739 735359 DAV2NPORT 5038 
= 51/7 7/22 2.3 tha BR2 SS 
= DIST PAINpH: 4.3 1067e8 
7 0 23 6.6 7.235 282 236 42 11 192 10 7 22 110 2.9 17 122 
2197 TUNIS LAKZ DELAWAKE NY 
421356 744649 ANDES 6088 


o1/f 7/22 1.08 Tha BR2 SS 
DIST3 RAINpH: 4.23 1092a8 


0 24 Jb2 336 294 251 45 15 191 29 7 26 129 1.45 24% 494 
2 23 8.6 9.03 296 231 53 14 200 113 6 25 129 1.7 23 528 
2198 SEXSSITH LAKE DELAWARE NY 
42275% 7449431 DAVENPOST 564a 
Ply TW/22 2.08 l3ha Ba2 SS 
XIST3 PAINoH: & 3 1)67a8 
9 %24 6.6 7.24% 129 87 37 12 130 113 8 24 61 47 «13 284 
3 24% wowed 7213 106 55 31 13 131 21 7 23 81 1.8 13 349 
2199 TRIB TO ELS CRESK ST LAWRENC NY 
442039 751256 SDOWARDS 225a 
bly 6/12 F8ten 2a BR3 S3 
"TST2Z RAINoOH: 4.2 101600 
0 4% 4.1 de52 006 606 79 23 305 26 25 31 251 9.0 5% 67 
297 LAKE GENEVISVE WAAREN NJ 
395937 745892 BLAIRSTOUN 1258 
sly v29 1.58 Bha BRS S3 
VISTS® RAINDH: 4, 3 121988 
9] 0 7.3 %F.2% 2349 2261 230 31 922 514 3 57 432 0.8 154 577 
1 4 72.4 $%9,.43 2571 2520 338 37 B46 543 4 77 454 0.9 135 S64 
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2210 SOUNTAIN LAKE WARRES a 
405130 74591) SASHINGTON 125a 
61/7 1/30 1.58 S2ha BR S5 
DISTS S3AINpH: 4&3 121988 
0 0 7.6 %e36 13937 3217 209 21 819 155 6 134 296 1.7 231 499 
1 &® 7.% 8.59 1757 1716 223 38 518 161 8 66 289 3.2% 127 509 


2211 LAK? SOLITUDE HUNTERDON NJ 
49940265 745315 HIGH BRIDG 87e 
617 1/730 2.38 Tha BB1 SS 
DISTS RAINpH: 423 121988 
0 09 7.5 68.39 1452 1383 228 16 625 389 4 169 183 0.8 300 344 
2 1 7.3 8.37 1507 3680 212 30 692 335 8 35 223 0.8 160 326 
z212 ALPKAUKEN CPELKR DAS HUNTERDON RJ 
4024622 745439 STOCKTON 508 
B1/ 1/31 1.28 6ha Ba2 SS 
DIST RAINpH: 4% 3 121988 
0 0 6.8 6.43 670 650 131 34 916 263 8 65 469 13.6 183 892 
3 3 6.3 6.99 629 589 99 350 4487 192 3 55 765 13.6 198 815 
oO 2213 HARIHOKAKE CREEK HUNTERDON NJ 
; a 403244 75094697 FRENCHTOWN 35a 
BIZ 1/31 e8¢98 se PR2 SS 


DIST@ RAINpd: 4&3 1219868 
0 9 7.2 62692 814 749 108 39 536 935 ) Le) 236 0.8 165 476 


221% 51% SIL? RON SO“SERSET RJ 
492812 743232 SONMOUTH JUNCT 1768 
31f </ll F00e Se Re2 S3 


DISTG PANDA: 43 121988 
0 0 6.9 6.92 3487 329 23) 52 244 639 10 a4 199 7.7 93 654 


2215 ROCK BUN POND SOSERSE&T NJ 
902499 7443153 BOCKY AILL 278 
uly 2f/i2 0.38 Sha B82 S4& 
DIST4 RATWot: 4% 3 121988 
0 0 6.6 7.2.1) 217 185 «6129 69 252 365 49 33 259 6.9 71 526 
2 0 6.4% 6.9% 231 199 116 79 383 %9 #849 33 449 6.8 61 511 
2216 PERRINEVILL? LAKE SONAOUTH NJ 
401333 74261) ROOSEVELT 508 
417 2/12 0.68 Sha 342 SS 
DISTS RAINDH: 43 1168088 
0 3 6.3 6.26 45 9 171 63 532 39) 1 a5 506 43 #77 236 
2217 R2D VALLEY LAKE AONACUTH NJ 
400946 7424835 ROOSEVELT 41a 
sly 2/712 see, Tha Ba2 SS 
DIST RATNpH: 4 3 168088 
1 % 6.3 Se3b 1771 lus 4? 84 162 ai3 3 39 229 12.8 $7 206 
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LITTLE PISHING CREEK CENTRE 
405322 773828 MINGOVILLE 


ol/ 3/3 
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LITTLE 


89/ 8/2 
DIST3 


1, ***2 
RAINpH: 
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SANDY CREEK 
412258 795543 POLK 


1 s¢¢g 
RAINDH: 
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LITTLE ALLEGHENY RES 
403110 782620 ALTOONA 


BR2 S5 
10 l6éaa 


VY ENABSGO 


BR2 S4 


1031a2 


BLAIR 


61/ 6/30 *¢e,— Sha BR3 S3 
DiST3 RAINpH: 4.2 10%36ma 
POWLEB HOLLOd@ ROW PERRY 
401558 773515 BLAIN 

380/ 8/26 ***s 38 BR2 S5 
DIST1 RAINpH: 4.3 1067a8 

N BR BOWSAN*S CREEK LUZERNE 
412123 761405 SWEET VALLEY 
50/ 9/23 #*ta 5a BR2 S5 
DIST! RAINpH: 42.3 965e0 
STRAIGHT RUN TIOGA 
414726 772417 ASAPH 

BO/ 7/24 *#8Qq 3a BR2 S5 
DIST2 RAINpDH: 4.2 965aa 
UPP2R THREE RUNS CLEARPIE 
411230 780785 DEVIL'S ELBOW 
olf VY 7 #eg 42 BR2 S3 
DIST2 RAINpH: 422 101608 
WYCKOPP RUN CAMERON 
411656 789800 DRIFTWOOD 

BIS T/16 =F#Fg 9a BR2 S3 
DIST2 RAINpH: 4.1 1) 16088 
RED AUN oLK 
411709 781439 DRIFTWOOD 

B1f B/17 eee 8a BR2 S3 
DIST2 RAINpDH: 4.1 1) 16m8 
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DISTZ PAINPH: 41 101600 
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2332 TROUT RUN LYCOMING PA 

412342 7703489 TROUT BON 2198 

B81f 7/29 eee 6e BR3 S3 

DIST2 RAINpH: 4.2 13 16m 

0 16 6.7 6.75 238 195 a4 6 166 1 3 16 71 #O.S5 8 195 

2333 SLATE &8UN LYCOGI NG PA 


413130 773151 LEE PIRE TOWER 3088 
B1f 7/29 ete 108 BR3 S3 
DIST2 KAINPA: 4.2 965a8 
0 116 6.9 68.33 333 259 53 13 201 34 2 24 123 0.6 21 143 


233% LEFT BR HYWER RUB CLINTON PA 
- 4912930 773602 SLATE BUB 3758 
a 81/ 7/29 seen 4a BR3 S2 
“3 PIST2 KAINDH: 4.2 965ee 
i 0 15 6&5 7.16 208 342 37 31 64 23 3 35 60 0.6 15 84 
“~~ 2335 COOK'S RUN CLINTON PA 
_~ 811819 775557 KEATISG 3358 
—_ 81/ 7/29 eee— 6a BR2 S2 
Pad DIST2 PAINpH: 4.2 96588 
fees, 09 7 %7.0 7.12 199 162 5) 23 117 7 8 32 99 9.5 22 145 
— > 2336 LONG PINE RON RESER ADASS PA 
N 395618 772615 CALEDONIA PASK 4158 
d1/ 7/31 6.38 S5Sha 883 S2 
DIST2 RAINDH: 4.3 101608 
0 25 S99 =%7.21 63 40 23 56 37 34 3 a9 47 61.9 16 y4 
1%6 23 S.3 5.83 18 -9 41 65 37 35 3 53 51 2.2 19 96 
2337 RED LION RESERVOL2 Your PA 
395633 703455 RED LION 1830 
oly 7/31 3.508 tha BR3 S2 
DIST2 RAINpH: 4.3 13 16a8 
9 26 7.3 7.25 532 Yel 1397 13 380 165 3 32 2856 1.0 66 601 
7 #5 6.8 5.89 1837 1805 191 15 745 205 5 26 462 eee 868 5 
233d OPOSSOA LAKF CUNBERLAND PA 


491432 771637 PLAIWPIELD 1378 
at7 7731 S38 2tha BR2 S2 
PISTS RAINDH: 4.2 10 l4eae8 


0 28 8.2 48.58 610 569 101 a5 257 119 3 45 299 2.9 9% 733 


7 1% 6.3 6.69 726 6437 120 33 323 114 5 37 326 «6888 665 890 
2349 LETTERSRENNY RESESVOI FRANKLIW PA 
#99051 774120 ROXBURY 2678 
31y 7731 & Ia 27ha B8L2 S3 
DiST2 CAINpH: 4.23 £Wil6ee 
09 2% 6.9 7.53 545 524 67 3 “89 4 5 12 169 0.6 13 96 
4 1% 6.9 7.23 609 572 78 3 99 33 6 12 157 3.9 13 143 
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FLD LAB FEP DEP 


DT pH pH ALK ALK COND Al Ca C1CLR K Mg Mn Na SO, 
23439 ScLAUGHLIN Ce2ee CRAWFORD PA 

$19018 7943233 TITUSVILLE WOR 4068 

RI Bf & Seeg 6a S8R4& S2Z 


DIST} PAINDH: & 1 1392m0 
0 18 6.8 8.55 1313 918 225 e 928 #1285 9 130-343 (005 «632 0—=«25 


23590 WOODCOCK CREEK CBAW FORD PA 
413953 7957548 TOSNVILL2 3030 
S1f/ df 4 8888 3a BB2 S2 


DIST RAINpH: 4.1 113967a8 
0 22 6.5 7.39 1165 i919 155 3 1075 192 4 1 330 22.1 27 288 


23517 TASARACK LAKE CRAWFORD PA 
413450 803442 COCHRASTON 3718 
b1/ Bf & 1.08 231iba B8R2 S2 
DIST3 RAINWpH: 321 94188 


9 26 #%F.6 8.93 425 369 97 4 s32 163 #11 4 177 #17 «7 174 
® 2% 6.3 7.31 5903 333 89 4 464 784 7 14 182 6.8 17 146 
2352 FRENCH CR PEK CRA@PORD PA 
819250 8098633 MSEADVILLE 3348 
. 81/ df & #0458 238 B8R2 S2 
DISTS RAINpH: 4.1 10S1e8 
9 25 %%.6 8.02 1429 1369 232 3 1696 316 10 1 529 3.9 42 a64 
. > 2352 PINE RUN S ERCER PA 
r > 411125 301155 sfecer 361.8 
sl/ s/f & 8008 ja BRS S2 
DIST3 RAINpPH: 42 Wiéee 
0 20 #=+75 %8&39 2682 2577 77 8 3879 632 5 18 #1792 1.8 235 2953 
15S GLADE RUN RESERVOIR BUTLER PA 
494257 79561) VALENCIA 1366 
31/7 37 4% 9.52 19ha BR? SS 
XLST2 S8AINpe: 8 1 965ae 
> %2£6 8.5 379 9239 865 173 5 1203 149 7 14 452 2.3 4%8 547 
* 20 6.83 %.19 #1606 1533 2191 5 1281 197 9 15 483 eee 590 281 
2355 PEWN BOOSPVELT LAKE CENTRE PA 
$04337 7749205 BARRVILLE s96ea 
Bis, b/ 5 S888 2s Ba2 S2 
DIST? RAINpi: 42 01688 
0 17 Se? 642) 54 18 15 103 25 5 3 68 #4 0.5 21 50 
2 33 Sed 6415 53 13 23 101 28 52 4 66 98 0.6 22 a1 
23457 8 BR STANDING STCNE HONTINGDON PA 
493855 776535 MPcALEVY'S FORT 2568 
sly 1/7 @ e089 53 BR? S2 
DIST3 F&AINpd: 4,2 965ea8 
0 14 6.3 7.16 95 63 35 Sa 61 38 3 a? 78 «6065 «626 153 
2358 GLOBE #0N HONTINGDON PA 
433653 775953 PINE SROVS SIL 27868 
417 77 #8 eee, 4a RR S3 
sTSTY RATNoM: 4.2 io5ae 
0 6 6.3 295 4h 8 ia 6 70 21 5 »4 96 29 37 127 








FLD LAB FEP DEP 


D T pH pH ALK ALK COND Al Ca cicLR k Mg Mn Na SO, 
2359 HIGH POINT LAKE SOSERSET PA 
394656 791350 SARKL2 TON 7568 
817 7715 .Je i tka BR2 S2 
DIST3 &AINpA: 423 T16faes 
0 25 69 $%7.79 789 707 366 9 559 57 5 7 102 1.8 5 778 
5 24 7.9 7.69 564 513 71 9 558 33 5 7 02 3.5 5 473 
2360 JORDON CREEK LENHIG# PA 
993717 753314 ALLISTON WEST 98a 
dif T/28 Shey 208 SSRs Ss 
DIST PAINpH: 4.3 1 18ea8 
0 23 8.1 3.30 3846 804 188 5 1279 332 7 15 “92 0.8 51 394 
<3©17 ANTIETAS SESER VOIR SPRAKS PA 
402124 755217 BISDSB8ORO 1598 
s 81/ 71/28 1.98 Tha BR Sa 
‘ DIST3 RAINpH: 7.2 1392a08 
—~> 0 2% %Y.09 7.59 4966 925 185 31 311 204% 10 36 389 0.9 39 10860 
= 10 19 6.6% 7.64% 962 929 156 38 344 207 #14 62 381 7.23 195 292 
we — 2362 CRABBY CREEL -HESTER PA 
4 — 400259 752803 VALLEY FORGE 98a 
- - B1/ 7/29 9088 3a BRS S2 
DIST3 RAINpH: 423 111888 
0 186 #+$%F3 «G27 1667 1898 368 3901805 366 5 36 06«6 12650 «60.9 (328 610 
2363 VALLEY CREEK CHESTER PA 
900301 753330 FALVERN 87s 
B's T/29 Ge4e6 is Bea ss 
DISTS PAINOAH: 72.3 318088 
0 17 #+%9.5 %%15 #18692 174s 3 23 2357 1529 2 310620710 «9.90 S22 yal 
26 PICKERING CRE2ZK CHESTER PA 
400695 753213 SALVERN 43a 
Aly 7/29 00a Se oR? ss 
DISTS RAINpH: 4.3 VIVAee 
% %22 %%5 8.43 4879 342 69 18 702 423 7 28 502 0.5 190 351 
2%%9 SARSE CRETK LAKE -HESTER PA 
4)0 343 754312 DOUNINGTOWSR 1108 
sly 7729 3,32 2Zivha BRi S3 
DIST) RAINpd: 43 992008 
0 26 8.9 47. 619 969 16) 12 552 609 5 23 333 0.6 %3 289 
9 20 6.7 7.37 S78 Sls ya 13 S74 ua2 5 24 3a5 11,7 #58 796 
2b YECA CASP PIND LESAWON PA 
HD 1445 Jo2gs7? SANNELS 26a 
siv 7729 1.58 tha 883 S? 
OTSTS PATSNods 406 1d67acn 
3 2s 7.3 4.94 582 545 46 92 125 9 11 65 atl 3-3 GF 1171 


Pa 22 6.9 4.99 665 630 97 98 92 115 13 »?7 176 1.4 83 101 
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FLD LAB FEP DEP 





D T pH pR ALK ALK COND AL Ca CL CLR K Mg Mn Na so, 
2376 SAW CREEK PIKE PA 
419758 7591355 TWELVEMSILE PON 3088 
B1/f 6/25 *ta 5a BR3 S2 
DIST2 RAINcH: 42.3 1168ma 
0 22 $.2 6.93 113 91 26 5 113 63 100 15 45 0.0 + 67 
2377 LITTLE BUSH KILL PIKE PA 
411023 750037 TWELVZSILE POW 3495a 
BIS 6/25 Fea Ja BR2 S2 
DIST3 RAINpH: 4.3 114308 
0 24 6.5 6,42 43 49 27 22 98 39 «660 32 739 «#721 «216 94 
2378 BRODHEAD CREEK MONROE PA 
410058 751155 = STROUDSBURG 1338 
81/ 6/25 **%g 6a BR3 S3 
DISTS RAINDH: 4 3 1219m8 
0 19 6.6 7.49 203 185 53 20 139 155 4 29 1171 0.5 22 122 
7a 2379 UNITY R®FSERVOIR WESTMORELA PA 
“£3 491650 79300) LATROBE 3378 
“4 B1/ 7/15 3.30 Sha BRY Su 
pee DOISTS RAINpH: 423 W)16e08 
eo 0 25 6.9 6.24 462 411 126 y 551 ih) ee) 21 299 0.0 41 900 
= 14 191 #+=%6.6 %7.08 462 407 131 3 588 wW5 3 20 299 15.6 41 965 
=P 
— 2380 KEDSTONE CREEK PAYETT2. PA 
= 395035 794222 BROWNFIELD 4 36a 
my me s1/ 7/15 F#eq 2a BR3 S2 
in DIST? RAINpDH: 423 109208 
0 16 6.6 6.71 227 166 58 29 143 432 3 34 130 9.4 32 1609 
2381 DEER VALL2Y LAKE SOMERSET PA 
394745 791135 MARKLETON 599a 
81/7 7/15 3.58 45Sha BR2 S2 
DIST3 RAINDH: 4, 3 liv3na 
0 25 6.4 5.32 105 51 43 3 13d 75 3 13 62 2.6 5 75 
5 23 6.3 6.291 139 7% 42 3 197 5) 4 12 60 2.4 5 163 
2352 NORTH FORK 22SERVOIR SOMERSET PA 
431603 790023 RACHELWOOD 4668 
Bi/ 7/15 4Ia 39ha BR2 S4 
DIST2 PAINDS: 4.2 114328 
9 25 7.5 3239) 524 483 74 21 170 #86 5 39 138 0.4 238 292 
5 16 6.4 6.73 168 125 63 28 149 57 5 35 135 0.9 33 314 
2333 CO@AN'S GAP LAKE POLTON PA 
499905 77553) BUPNT CABINS 3718 
BIY 7/11 2.8m 290ha Ba) S2 
DIST3 KAINDH: 4&, 3 9+ 5aa 
9 27 64.2 4.77 #4219) 67 33 42 71 53 5 44 80 10 24 ¥6 
4 17 6.9 6226 19) 67 34 uS 67 53 8 44 79 6.4 25 96 
al 
2344 oP AVEB SEL! FSSERVOIP CA@APIA PA 
4¥O%3629 784191 BEAVERDALE 728a 
RIS T/ 6 1, >a edna RKS S2 
QISTVT RAINS 4.2 111#ea 
6 22 S7 6.33 31 4 27 29 54 21 3 35 51 2.4% 12 5? 
3. 5.5 6.45 45 13 27 32 Ib 23 3 37 53 2.5 18 46 











FLD LAB FEP DEP 








D T pH pR ALK ALK COND Al Ca CL CLR K Mg Mn Na SO, 
2385 SPORTSSIAN'S LAKE BEDFORD PA 
395758 7829357 CLEARVILLE $37a 
B1/ 77 38 2.58 6bha BRY SS 
DISTS RAINpH: 4, 3 965aa 
0 27 #69 #7.59 158% 22 a9 33 156 39 3 37 158 0.9 42 163 
6 16 6.9 7.98 268 220 73 27 187 65 3 35 172 3.2 43 304 
36 LAURZL HUN HONTINGDON PA 
4$02015 783633 CASSViLLE 824e8 
81/f lf 8B 8#t_ 32 BR2 S2 
IST2 @AINDH: 422 W)16a8 
0 18 5.6 5.42 53 © 348 266 27 33 2 Wo 76 #0O.5 61 1648 
= 2337 SONPISH POND BRADFORD PA 
4133460 764148 LEROY 632s 
81/7 7/7 7 2.38 Wha BR2 S2 
DIST2 RAIWpH: 4.2 I1%e8 
0 25 =%7.4 6.16 49 9 2 4 118 35 3 ta 46 2.3 4 33 
3 23 7% 5.95 4? 22 24% 4 126 31 3 13 au 2.5 4 78 
2388 OLYa2TC LAKE SUSQUEHANNW PA 
415559 7o57S18 FRIENDSVILLE 5038 
“= 31/ 7/7 2.98 Wha BR3 S3 
= DIST3 RAINOH: 42 918a8 
p= 9 25 7.1 7.41 172 185 393 9 142 42 4 21 76 64.7 ~=«10 101 
no 5 22 7.2 6.93 149 193 36 19 166 ee 93 47 14 101 
2389 TOSCASOKA LAKE SUSQUEHANN PA 
414420 76952) AUBURN TENT2® 3538 
sly 7/7 7 9.58 25ha B23 S3 
DIST2 FAINDE: 422 918e8 
59 25 10.3 3.393 299 262 55 11 129 63 17 22 71 18 #17 94 
5 13 6.5 6.36 592 u71 79 9 296 u3 8 20 193 16.9 22 35 
23490 LAK® J®Aa LUZERNE PA 
4120909 76180) R2ZVaRICK 478e8 
31/7 7T/ 7 3.33 WWOIha BR2 S2 
DJISTS BAINOH: 4.2 19 l6ea08 
09 22 %5 %.21 2? p) ‘) 0 61 23 3 rs 13 «(1.8 0 89 
4 23 5.5 6.66 31 2 29 9 W5 52 3 5 yw 2.3 0 122 
234917 SILVER LAKE PIKE PA 
419152) 745845 EDS ERERE 3998 
dif 6/25 5% 58 S7Tha 884% S2 
DIST? RAIWoH: 4%. 3 111888 
0 23 6.1 6.42 49 13 26 8 11 26 + 1 S55 1.0 7 111 
5 23 6.8% 6.57 u5 18 2b 9 105 3 J 29 53 1.90 A 101 
2392 DES? LAKZ NONaU® pa 
4103190 39592359 POCONO PINES S5/ea 
81/7 6/25 3.58 Yha aS S2 
OIsTy RAINpDH: 4 3 134'sea 
J 22 962 4,35 22 J 23 6 54 34 1 16 28 3.9 ; Ju 
ry 24 Sed 4,35 3 9 <3 6 61 iu j 17 27? 3.9 } 62 








FLD LAB FEP DEP 





D T pH pR ALK alk COND AL Ca ch CLR EF S& HE RS 80, 
2393 SHADIG®?® CREEK e AYRE PA 
915353 752783 STABBOUCCA 4098 
b1/ e/22 eeeea 138 BR3 S2 
DISTS KAINpDH: 4.2 106788 
0 20 #%%7.8 %73.3/5 434 398 83 0 aa5 182 #610 38 95 1.3 5 89 
239% LONG POND @AYSE PA 
4143000 752037 ALDENSVILL® $138 
31/ 6/22 43.58 33na B8R2 S2 
DISTS BAINDH: 43 tIl11éae 
0 2% #%79 +744 199 172 89 9 274 52 10 9 64 0.9 4 132 
0 2% =%+F.9 6.98 29% 267 69 i) 121 73 2«s«:119 7 24 5.8 1 101 
2395 @ BR DYRZRBY CREEK eAYNE PA 
413937 751725 ALDENVILLE 3178 
31f 6/22 e888 10a Be2 S2 
— MIST2 RALSpH: 4.3 1119888 
> 21 #=+%F9 %$38.655 483 853 89 0 360 138 «610 6 77 #+O.7 4 117 
+ pao 2396 DOYSERRY Cee WAYHE PA 
4 > 413583 751632 HONESDALE 2968 
~ ~ 51/7 6/22 eee, 158 Ba2 Sz 
— DiSTS BAINpH: 4.3 1180828 
wo 0 241 #=~F9 #$%$93 #443 311 81 n) 892 121 10 8 99 1.1 5 137 
- 2397 LACKAWAXEN SIVER WAYNE PA 
- 4134335 751615 HONZSDAL® 31628 
B17 6/22 e883 258 Be2 S2 
DISTS PAINoH: 4 3 1131888 
9 23 7.8 8.392 348 312 77 ») 397 142 + M 98 2.6 ) 153 
23396 LAK® LACA@AS WAIWE PA 
412236 751731 HAWLE 4398 
B1/ 6/22 3,58 2tha 8hk2 S2 
DISTI RATSpH#: 43 116888 
0 2% 6.9 6.04 9 77 35 2 168 ‘4 10 11 49 6146 3 106 
114 120 «548 «5S «6100 67 5% 0 73 23 19 9 44 23 2 84 
2399 @ALLINPAUJPACYK CRESEK PIKZ PA 
411355 751855 NEWPOUNDLAND 3938 
41/ 6/22 S889 158 RR? S2 
DISTs PATNoY: & 3 1163aa 
0 19 7.6 7.9%) «6285 249 67 3 250 134 5 19 79 «#2609.5 6 101 
4491 CUBBAGZ POND SOSSEX pF 
1°5156 752339 ELLENDALE Qn 
agysvity 7 1.98 tha BR’ Ss 
OiSTI PAINpH: 4.5 119888 
0 11 6.5 be32 537 104% 13 2b 239 271 13 4 360 COW O88 57 
2 10 6.2 7.205 342 t14 %5 33 249 253 iad 5d 115 9.8 159 52 


103, 








FLD LAB FEP DEP 


D T pH pR_ ALK ALK COND AL Ca Cl CLR K Mg Mn Na So, 
2402 SILVER LAKE 2 NEW CASTLE DE 
393726 753332 SAINT S2ORGES Se 
a0/vyilly 4 90.58 lgha Ba2 Ss 
DIST3 RAINDH: 43 111888 
0 13 7.9 7.37 832 39% 192 66 285 34) 5 32 285 5.2 73 309 
3 1% 7.0 6.3% 43935 389 «6123 76 302 337 12 33 338 5.17 648 314 


2803 LUAS POND NEW CASTLE DE 
393345 754357 SAINT GEORGES 13a 
sosiiy & 0.58 77ha 52 S4& 
DIST3 FAINpH: 4.3 111888 





0 10 6.9 7.26 480 $32 13) 41 632 385 2 62 389 1.7 117 452 
1 11 6.2 7.47 489 437 1285 35 651 343 3 39 345 1.2 141 432 
2404 SHALLCROSS LAKE NEWCASTLE DE 
393054 754925 MIDDL&TOMN 12a 
= B9/1t7 & O.5a 19ha BR2 S3 
mr DISTS RAINDH: 4.3 111888 
S 9 10 7.0 7.19 722 637 176 5S 839 #343 8 42 661 4.3 90 361 
_ 3 10 6.6 7.09 632 579 195 61 1073 353 5 39-973) «3 B2 364 
=> 
= 2405 SORRIS MILLPOND SUSS2x DE 
as 383844 751654 HARBESON 7s 


aoyvyilly 7 9.98 Sha 8RB3 S3 
DIST3 RAINpDH: 45 111888 


UTI 
Yo) 


6.9 6.59 332 308 91 45 121 311 2 Su 76 2.6 199 162 
a 1 9 6.2 6.738 361 328 89 41 165 396 «25 62 101 2.6 120 1490 


2406 SILVER LAKE 1 NPW CASTLE DE 
392616 754222 AIDDLETOWS 12a 
dJsiiy & 9.58 tha BR2 S3 
DIST3 PAINDH: 42.4% 1113388 
7 72.23 006 579 211 49 585 538 6 a5 424 10.3 99 614 
3 19 65 7.36 613 566 211 53 593 551 4 49 468 12.1 93 611 


2407 HAVEN LAK® KENT DE 
3354866 752733 NILFORD Tt 
#0/1l7 6 O.Y8 J3ha Bad SH 
DIST® RAINpH: 45 111888 
0 11 #63 6.97 366 337 99 36 407 3350 20 58 215 3.4 136 202 
12 6.7 6.91 4323 295 192 54 349 311 20 u1 275 3.4% 92 216 


La) 


2408 TRAP POND SUSSEX oe 
383126 752823 TRAP PUND 1200 
80/117 6 1.178 S3ha BRI S3 
DIST} RAINpH: 4.5 111808 
10 6.4 6.33 19) 12 5S? 45 34 «61938 5) 23 


J «7 199 119 
1 10 6.2 Se62 19) 157 57 41 73 213 #10 5 45 7 


117 145 


—_— —_ 
* 


2399 LAKE COMO K ENT DE 
$¥1730 753613 SSYRNA Sa 
scsvily S&S 3.38 Wha BR SS 
DIiSss WAINpDH: 4.4 11 18es 
4 51 439 
0 55 473 


) 9 Se9 6.651 313 285 114) 97 353 335 8 23 506 10. 
2 WW 6.4% %.96 342 309 1348 38 364 343 1 27 446 11 


10 








FLD LAB FEP DEP 
ALK ALK COND Al Ca cl 


io 
i“ 
~~ 
x 
fm 
QO 
re 
mm 
ix 
be 
x 
3 
- 
a 
“4 
° 





2410 GARRISONS LAKE K EST DE 
39 1523 753522 DOVES 3a 
s0/livy 5 3.58 35ha BRI SS 
DISTS RAINpH: 4.4 1113888 
0 9 7.0 7.19 731 660 174 90 587 377 «13 30 793 6.1 58 350 
1 9 7.2 7.19 674% 641 to7 71 586 417 618 42 601 6.8 69 458 


23117 MASSEY'S SAILLPOND KENT DE 
391821 753732 DOVER 5a 
B0/11/7 5S 0.58 Sha Ba SS 
DIS14S RAINpH: 4. IT11688 


9 10 7.9 7.37 753 718 #224 Se 666 383 10 49 540 5.1 86 $23 
1 10 7.0 7.13% 741 £4793 23% #=+%SB 773 4427 6 42 671 6.0 6&8 523 
2412 CAIPSAN POND SUSSEX DE 
383345 753220 LAUR2L 6a 
: 80/117 6 1,18 13ha BR3 S3 
a DIST3 RAINpH: 4.5 111888 
> 0 110 6.4 6.95 260 252 94 31 418 274 15 74 199 0.8 158 202 
1 19 #67 #+%7.19 285 #257 #42x96 #$=$%29 +1382 277 18 28 163 1.7 168 223 
“a — 2413 FLEETWOOD POND SUSSEX DE 
= 7 383838 753032 SZAPORD EAST 9a 
= BO/11/ 6 See— 128 BRI S3 
» DIST3 RAINpDH: 3.5 Iliidee 
— 09 410 6.1 6.39 276 248 #61 32 95 213 «65 87 55 1.7 186 93 
x 
— 2414 SUD AILLPOND K ENT DE 
390613 7544656 MAR YO2ZL 1Je 
BOsI1/ 5 ete6 Sa = BR3 S5 


DIST3 RAINpH: 4.5 111688 
0 12 68 6.38 319% 295 128 49 312 300) 3 8 225 147 792 419 


2415 WYOMING LAKE KENT DE 
399740 753320 DOVER 9a 
09/117 5 96.58 2tha BR S4& 
DISTS RAINpH: 4% 1118688 


0 10 6.9 6.245 379 3370 = 481 $2 314 361 8610 62 193 2.6 118 335 
1 10 7.2 6.94% 379 342 152 42 399 377—S 45 54 250 3.4% 115 348 
2416 VOSHZLL SILLPUND KENT DS 
390535 753146 WYOMING 6e 
BO7117, 5S O.5e 1S5ha Bas SS 
DIST3 RAINoH: 4&4 11 l#e8 
3 11 6.9 6.5) 333 318 35 au 432 427 + 59 257 4.2 112 u3u 
2 10 7.2 6.55 £323 290 ley 43 3%) 344 9 58 259 4.3 190 416 
2417 ANDREWS LAKE ¥ ENT DE 
3990130 753032 WYOMING Ys 
HO711/ 6 9.68 Jna S8a3 Se 
DISTS RAINpN: 4924 £111868 
fv) 9 6.1 6.47 276 257 3 41 246 365 10 53 149 6.9 11? 369 
; 9 6.1 o.77 297 295 #937 a4 2389 39> 14 52 185 6.0 110 369 











FLD LAB FEP DEP 
D T pH opr ALK ALK COND Al Ca cl Cc 


tr 
~ 
= 
be 
x 
td 
f 
e 
vw 
°o 


2418 MARSHY HOP® CREEK KENT DE 
385058 754025 SICKSAS 128 
B0/11/ 6 eee 7s Be3 Sse 
DIST3 RAINpH: 4.5 111888 
0 8 6.4% 6.32 261 228 38 a9 365 279 3 56 266 5.2 99 366 


2819 KILLENS POND KEWT DE 
385753 7531546 HARSINGTON Sa 
so/117 &© O.o8 JJna BE? SS& 
DIST3 RAINpDH: 4.3 1118288 
0 9 6.5 6.39 285 257 153 a5 a79 372. 2 55 324 S22 197 S93 
2 6.6 6.59 291 266 159 32 692 377 _ 22 62 39 43 113 506 


ar 2420 COLLINS POND SUSSEX £ 
334236 753123 SEAPORD EAST 9a 
so/11/7 © 9.98 Jdha BRI SS 
DIST} PAINpH: 3.5 11718m0 





9 #+$+.9 %5S.9 6.32 181 147 51 ua 77 lo 8 55 50 1.7 110 149 
2 9 6.0 6.48 19) 162 58 50 113 185 3 a4 82 1.7 $9 166 
-~ 2927 HEABNS POND SUSSEX DE 
ho 389041 752534 SBAPORD EAST Sa 
vosily 6 O88 2Zjha 8R3 SG 
DISTS RAINpH: 4.5 1319888 
- 10 6.3 5.665 37) 336 4 T18 35 307 3938 4 68 206 2.6 126 299 
1 10 6.2 6.71 366 326 #111 39 369 298 i] 81 187 0 63.4% 139 333 
2500 2 BLAIR VALLEY LA*2 SASHINGTON 8D 
394156 775635 CLEAR SPRIWG 2108 
Bly 2/17 1.58 Tiha BR2 SS 
DIST} RATNoO4: 4& 7 965aa 
0 9 6.7 %27.5) 44% au07 89 is 557 63 21 67 263 0.0 1448 390 
8 S 6.9 7.19 70? 680 136 37 779 59 5 76 #20 90.0 131 %32 
2503 ROCKY GA RUN ALLEGANY aD 
#494257 743336 EVITTS CREEK 57a 
HOsSI2Z/27 «Fete 23 BRa $3 
NIST2 FAINDH: 4, 3 914e08 
0 0 6.7 6.61 13569 1406 18) 34 184 147 3 74 94 0.0 141 385 
2509 PINEY CRE -K sARRLTT ad 
4942186 7999455 AVILTOS 7118 
AVs/12Z/27 Stee 43 PRi Se 
DIST2Z RAINpH: 43 116888 
3 1 6.2 % 83 £4185 I4e «6107 55 37 42] 4 44 350 1.7 #49 241 
25905 CRABTRES CREZK SAPSRETT Av 
144916 79092) BITTINGER 4608 
» 412727 eee, 53 Sez Se 
OIST2 waAingu: 4.3 121988 
0 1 bel 6.17 435 jas 124 4? 634 lds 5 41 #29 9.9 1935 377 


M2 





FLD LAS FEP DEP 


D T pR pR alk alk COND AL Ca ch ciR *£ “gs we SS 30, 
2506 LAKE SABEES ALL EGASY aD 
394205 763912 EVITTS cCREer 3560 
gif 2717 1.58 79nha Bes Se 
DIST2 RAISp#: 4 3 3 14a8 
i] 0 7.3 8.1) 615 S76 83 &6 625 35 3 85 552 1.7 196 339 
15 '* 37.2 8.32 $53 518 95 39 633 e) | & 92 283 0.9 162 343 
2507 FROSTSI&S RES®RVOIR sAR BRETT aD 
394249 790020 AVWVILTOS 7178 
817 2717 2.38 i2ha BRI S3 
DIST2 PAINp4: 4 3 ilese8 
2 09 6.3 6.45 126 62 6? 3a 1é6 219 4 74 98 9.9 1 274 
5 2 6.3 6.35 186 Je «4458 55 337 ee80 5 as 350 1.7 89 351 
2508 ShAVELLY PUN SARRETT aD 
3932198 792082 ScHeNRY 7S3s 
BI/IZs/27 fete 2s BR2 S3 
D'STI RAINoe: & 3 127088 
9 0 6.1 6.8) 196 93 89 27 139 35 i) 92 55 0.9 180 208 
2509 GOSPOWDER PALLS BALTI®ORE apd 
398126 764650 LIN ESSORO 183e 


BI/I2Z/29 Fete de Ba2 SS 
DIST2 RALWoH: S&S. 3 11888 
i) 0 6.8 ¢€.92 $27 $81 108 3? 598 39 3 70 336 0.9 139 105 


Lt 


2530 488aTSAN SSANCH BALTISORE AD 
- 392738 763521 TOWSON 9168 
B0/12/29 e886 2a Ske SS 


DIST RAIWoH: 43 116888 
0 1 69 %%.13 7% 777 Wa 31 954 768 3 99 $69 0.0 156 239 


251317 Calas ey BRAWCH BALTLAORE AD 
392326 76513) RFISTEPSTOEN ite 
BO/12729 eee 2a BR2 SS 
DISTs RATSpH#: 6 3 11138808 
*) 1 %7.5 7659 28612 2769 238 20 21 a5 j 27 37 9.0 239 216 
2572 NORTH PORK CREEK PREDERICE a) 
592746 7713121 LIBSRTYT OWN litle 
BOsI12/29 Fees oT Bee Ss 
DISTS® &AIN08: 43 1067e08 
0 1 %.0 7.28 21 334 dt 16 a7 232 0 71 262 9.0 136 93 
2512 OWENS CREE PREDZRICE AD 
19356 772008 wWOODTBORO ade 
S0/12729 886 5a BRI S4 
MOiSTS HALNpMs 4. 3 WOSle8 
0 0 %Fe% $%Foe98 13419 BVT208 297 40 995 345 3 93 44 08.9 140 S13 
257% HUNTING CHEER PEE DERICK SD 
193752 7727786 BLUE RIDGE SUR 2438 
49712729 ees ja BR. oe 
DISTZ PATWNpH: & 3 1) lhee 


0 5 7.9 72.3% SOT #54 )4 149 “51 216 3 57 257 0.9 128 218 











YLL 


2515 


2516 


2517 


2518 


2519 


2529 


2522 


"1 
tv 


D 
PATUXENT RIVER HOWARD aD 
3393417 770323 SANDY SPRING 11i2a 
B9/12/30 se*, 12a Ba2 SS 
DIST RAINpDH: 43 1067an 

0 
HAWKINS B2IVER MONTGOMERY AD 
391026 770120 SANDY SPRING Qua 
o0/12/30 tty 72 BR2 SS 
DISTSY RAINpH: 4.3 067ma 

0 
SILBERT CREPK CHARLES aD 
3830945 7638528 HUGHESVILLE 33a 
BY/12/30 eeey 3a BR2 $3 
DIST3 RAINpH: 44 13 16aa 

0 
TEINITY CHURCH RUN CHARLES MD 
332719 765199 CHARLOTTE HALL 15a 
81/12/39 eta 3n BRY S3 
DIST3 RAINDH: 4.4 1) 16m 

0 
TRIB TO OCTORABO CRE CECIL MD 
394119 76682) CONO#INGO DAS 33a 
BIf 6/7 1 Fete 3a BR2 S5 
DIST SAINpH: 423 1118a8 

0 
LAKE FOSHARP GARRETT AD 
392705 792820 OAKLAND 731a 
olY 2/17 1.56 Sha Ba3 S3 
DIST3 RAINpH: 43 132108 

0 

2 
PLUS CR2EK CECIL D 
39 3405 755310 NORTH BAST Ja 
d1f bf 1 F¥#~a 52 BR2 S3 
DIST3 RAiNpd: 4.3 111800 

v 
CYPRZSS BRANCH KENT 4D 
391525 75590390 MILLINGTON 32 
81/f 6/7 1 &*#~e 5a BRY S3 
DIST3 PAINDH: 4&4 1llsoe 

) 
@iLLISTON LAKS CAROLINE AD 
344940 75503) HOBBS 38 
s1f 6/ 1 1.22 l3ha Bh2 SI 
DISTS® RAINpH: 4.5 11808 

9 

3 


nu Oo 


17 


22 


23 
20 


5.8 


ww 
o . 
Wir 


6.7 


nn 
« 
oO 


LAB 
pr 


6.59 


6.04 


6.61 


3.900 


6.57 
6. 35 


7.53 
6.61 


FEP 


ALK 


263 


370 


194 


3 33 


6 8z 


31 
624 


2) 


353 


304 
313 


DEP 
ALK 


227 


336 


301 


848 


596 





76 


96 


98 


¥9 


148 


38 
35 


42 


88 


ohe) 
91 


39 


39 


45 


53 


49 
71 


993 


115 


301 


473 


472 


262 


188 
218 


38 


251 
163 


2u8 


319 


211 


163 


266 


31 
53 


100 


33) 


w 


70 


40 
37 


\* 


7i 


68 


60 


46 


75 


51 
ete% 


65 


72 


35 
51 


be 


197 


172 


320 


364 


558 


136 
232 


66 


0.0 


1.4 


5.8 


—_~ — 
6 
Nw 


124 


126 


199 


91 


296 


100 
63 


30 


55 
89 


so 


61 


182 


464 


408 


300 


232 
271 


189 


199 


179 
137 








FLD LAB FEP DEP ; 
D Tt pH pR ALK alk COND AL Ca ct CLR & ss Se %F so, 
252% HIGSINS S’LLPOND DOFPCHESTER &D 
3683119 7557465 TAST NEW SARKET 20 
81/7 6/7 1 o.3s 63ha BR2 SS 
DISTS RAINpH: 4.5 116808 
0 24% %F.2 6.99 365 416 12) 78 210 417 100 59 319 41.7 129 268 
2 2% 7.1 6.92 4481 397 #117 75 213 359 100 57 we 3.4% 123 205 
2525 aDKINS POND #ICOSICU aD 
3382999 75223) BANGS 9a 
ots 6/7 1 0.>a Jha 882 S48 
DISTS RAINoH: 4.5 121982 
0 29 6.6 6.42 333 296 196 73 134 329 79 55 1971 1.4 #71 242 
Pi 17 6.3 233 441 3392 194% 33 189 3430S 60 38 188 131.7 4&8 1862 
2032? SR2CH2ZS BIN SARSSALL wV¥ 
$35335 AQVI3II8>D FOUNDSVILLI 2358 
gif 3/719 tts 58 Bei se 
DISTS RAINpi: 4.2 1967eaa 
) 2 7.2 79.93 1377 1338 25) 0 1969 132 3 lv 515 2.35 36 927 
t 2003 CONAWAY LAKE TYLEE hi 
392520 895139 SHIRLEY 2j\e 
HiY 3719 1.43 Tha 883 S% 
DISTS S8AINDY: 3, 2 ll iseas 
= 0 2 7.9 7.28 441 492 19? 20 314 37 6 293 246 2.3 52 452 
- 8 3 6.8 5.95 431 997 83 18 352 63 5 26 260 2.5 51 459 
=: 2610 SAST LYNN LaKt WAY av 
160490 S22300 BAYN? 1968 
7 81/7 3/29 3.78 SIIha BRI S3 
— DISTS RAINpH: 4.3 111608 
ww 0 % 6.5 7.27 142 116 57 259 71 47 3 1d6 194 3.4 169 #95 
3 4% 6.7 7.99 137 193 69 261 71 73 4 135 192 4.0 185 a17 
2612 HORSE CREEK LAKE @YOSING eV 
373425 3814235 RAILOYSVILLE 3552 
stv 3720 0.528 Sha BR2 S3 
VIiSTZ EAINDH: 4.9 1118es8 
0 3 6.7 6.95 191 147 16 162 74 16 1 “9 1786 2.6 111 935 
6 3 6.6 7.79 296 166 76 44 43 25 9 a. os 614C 3 O15 439 
6135 PLUS ORCHAS)D LAK FAY2TTS wV 
375700 619439 SOAK HILL 5 28a 
JI4 J21 4.22 Iiha BR2 S3 
NISTS PAINOY: 4, G 1163288 
% 6.8 7.9009 167 138 45 138 77 2) 3 77 154 145 82 318 
8 3 6.8 6.79 #4167 143 51 199 33 36 3 68 1990 1.7 75 329 
261% BFONCOVE LAKE SONROE av 
5737)3 BV2107 PAINT GANK /63a 
S17 3771 1,59 She BRGY S3 
;3 O°YINeH: 4. 4 lJ4laa 
9 2 6.7 h#.51 147) 113 33 37 4 »5 » iy I 1. 3 21 127 
2 €.8 6.54 182 198 10 17 i) 45 5S to 60 1%? 24 127 
(1S 





FLD LAB  FEP DEP 


DT eH pH ALK ALK COND al Ca clcLR xk Mg Mn Na so 
7615 SHERWOOD LAKE GREENBRIER WY — —és 
380610 89903) LAKE SHERWOOD 6132 
olf 3/21 1.02 69ha BR2 S3 
DIST3 RAINDH: 4.4 1) %6e8 
0 3 6.8 6.225 73 39 23 52 36 35 3 a6 96 3.17 15 67 
5 3 6.8 6.72 193 73 20 0 65 29 6 10 17 ty 1 65 
2616 SPRUCE KNOB LAKE RANDOLPA wv . 
364210 793515 SPRUCE KNOB 11712 
Biy 7711 1.08 liha 882 S3 
DIST3 RAINDH: 4.3 Wtl6aea 
0 25 %97.5 7.34% 165 143 27 4 154 7 10 13 54 2.3 5 41 
5 16 6.5 6.4% 491 356 a5 9 133 5 € 2u 72 29.5 190 15 
2619 T2TER CREEK LAKE BARBOUR wv 
3399625 795239 YONTROSE 5608 
81/3/22 1.5m 13ha BRI S3 
a DIST2 RAINpH: 43 132188 
0 3 6.6 6.73 6) 138 35 78 53 19 1 56 82 2.3 32 156 
5 3 6.6 6.58 65 19 31 66 48 25 1 54 69 1.9 26 182 
— 2621 LITTLE RACCOON CREEK PEESTON wv 
- 392235 795010 NEW30RG o298 
a 41/ 3/22 ea 4a, «6B S3 


DISTY RAINpHd: 31% 127088 
0 3 Fed %Fe2db6 255 296 59 7 274 253 3 17 125 %3 #15 333 


TIA y 


2622 LAKE Of WOODS PRESTON we 
394295 793355 BRUCETON BILLS 571la 
BI1/ 3/22 SF%e 29ha BR? S3 
DISTS RAINpH: 4? 121988 





0 2 Fe1 6637 273 250 71 9 396 41 ") 9 105 1.5 7 318 
2623 WN BR SNOWY CREEK PRESTON wy 
392715 793109 TERRA ALTA 777a 
HI/Y 3/22 e483 6a BR3 S3 
DIST FAINDH: 4.3 132188 
0 3 7.9 6.30 uj 63 5) 17 690 au 9 27 59 0.9 11 261 
2624 STUNZY BIVER LAKE GRANT a’ 
399412 791755 MOUNT STORS 19 358 
dis 7711 3.72 %65ha BRI SS 
DIST RATINpH: 43 118808 
9 £3 Sed SI) 13 -y 22 19 47 13 3 28 35 2.6 6 67 
9 18 5.7 4&, 38 19 -9 23 23 44 10 5 32 35 2.8 7 84 
2627 CACAPON PARK LAKE SORGAN wv 
393090) 741895 RIDGE 3118 
s1/ T/11 2.68 2ha B3R2 S5 
JIST3 etAINpHA: 4, 3 91588 
0 26 7.2 6.93 977 947 116 3 946 34 3 12 313 3.9 27 117 
5S 15 6.9 7.31 1274 $%$+w140 1371 0 1926 37 5 11 285 36.6 21 72 


Ne. 











Lil 


2026 


2630 


«633 


2634 


Lae 

Co 
a 

Ww 


2636 


203? 


~0e9 


> tI 
WARDEN LAKE HARDY av 
399755 783550 YELLOW SPRINGS 3928 
61/7 7/11 3.308 18ha 8R3 SS 
DISTI PAINpd: 4%. 3 918ea 
0 27 
6 18 
SIGNEE RUN RITCHIE uv 
399455 81994) SMITHVILLE 2108 
B1/ 3/19 See ‘ea BE3 S5 
DIST3 RAINpaH: S22 1219288 
a) 2 
RT FR BOFFPALO CREEK LOGAN av 
375335 629949 CHAPMANSVILLE 2138 
B1/ 3/20 eee, 3a BR2 S3 
DIST3 RAINpH: 4.2% 1118m8 
0 4 
STEPHENS LAK2 RALEIGH wv 
374649 81183) ECCLES 6178 
81/ 3/20 4.00 93ha BR2 S3 
DIST3 RAINpa: 3.4 116808 
4) 3 
12 3 
PAUTHER CREEK NICHOLAS wi 
381325 803948) NETTIE 6778 
B1/ 3/21 eee, 38 BR2 S3 
DiST3 RAINpH: 4.4% 121988 
0 2 
BIG DITCH LAKE WEBSTER eV 
332410 803400 COWEN 6778 
817 3/21 L390 us BR2 S3 
DISTS RAINpH: 4.4 124408 
0 1 
5 
PICKLES PORK RESERV BRAXTOW WY 
384550 303555 ORLANDO 259a 
b1/ 3/22 *e¢e Sha BRI S3 
DIST3 FAINpH: 424 121988 
0 1 
STONECOAL LAKE LEVIS av 
365645 832250 ROANOKE 34a 
81/7 3/22 1.Jm 239ha BRI SS 
DISTS RAINpH: 4.3 121908 
") 3 
3 3 
IPPER GLADE C2 R2SER FALZiGH aV 
3743020 419345 SMYADY SPRING i238 
BIS 4/25 808% 33ha BRe S3 
DISTS RAINDU: 4.4 1119808 
> 12 


6.7 


6.7 


ae 
ee * 
~ & 


~~ oO 
*- * 
oc 


LAB 
pa 


60. 43 


7.02 


FEP 
ALK 


670 
1688 


465 


134% 


109 
137 


274 


142 
166 


490 
534% 


DEP 
ALK 


606 
1615 


832 


76 


78 
107 


235 


108 
118 


515 


456 
uy¥5 


}3 


94 
185 


112 


136 


33 
32 


239 


43 
111 


104 


113 
114 


14 


120 


64 
36 


12 
11 


56 


716 
1553 


202 


257 


61 
62 


439 


109 
11 


290 


236 
263 


70 


18 
15 


166 


24 


34 
39 


1433 


125 
625 


47 


167 
158 


612 


@w 


\* 


18 
11 


25 


74 


93 
39 


13 


23 
22 


20 


ls 
ly 


Sv 


334 
435 


126 


498 


37 
68 


155 


65 
64 


155 


104 
139 


91 


48.5 


0.5 


2s 
* e 
ww 


oC oO 
*-* 
2 © 


42 
33 


22 


246 


32 
18 


13 


21 


7 
19 


186 
67 


429 


856 


212 
202 


137 
94 


377 


$23 
516 


297 





FLD LAB FEP DEP 
D T pH pH ALK ALK COND Al Ca cicLR k Me Mn Na SO, 
2701 PEACEFUL VALLEY LAKE FREDRICK VA 
3908460 78621485 WAYPIZLO 2598 
oly 872% 1,58 Jha 3882 S3 
DISTI RAIWoH: 4, 3 3 ises 
3 16 7.9 3.63 975 936 111 0 1439 3) 3 j 209 9.9 3 228 
s 15 7.7 3.36 1146 11317 #47 0 1304 a4 5 295 3.23 16 273 
2702 SILL CREE SMSNABDOAH VA 
3445290 7564915 COSIZCOVILL®E ~Slea 
SIs 4724 F988 7a BR2 $3 _— 
- DISTS RAINpH: 6.3 914a0 ‘id 
0 6 8.1 %2F,38 199% W989 282 16 2625 199 5 26 1691 O.9 324 243 
2703 DEFY RON ROCK INGHAMS VA 
36246839 76833123 ELKTON BAST 6158 
GIF 4/24 8a 48 SRS S3 
— DISTS RAIWpH: 4S 1) l6a8 
= 9 13 #73 8.39 38? 356 52 s7 34 25 3 ud 111 0.0 38 52 
S 2793 LAK? ALSESABLE ALS IARL® VA 
— 309525 78373) CIOZET 152s 
rv sly #723 1,58 Wha #82 $3 
-_ DISTS FAINpd: 4.4% 91688 
-_ 0 18 72.5 6.43 338 313 SJ 16 26 106 ¢ 2d Sb 9.2 146 52 
x 6 0 6.7 7.98 338 289 5 25 195 39 % 32 39 6463 = «220 57 
<795 FLUVANYVA SURI TAN LAK PF LUVANWA VA 
375325 Fee2Is BOYD TAavVeRN lite 
Bly as24 i. J3 23ha BR2 SI 
oISsST3 RATA@OH: «4 5 llédee 
9 9 %F.%4 4.63 £52 224 3308 43 52 6 1 195 9.0 63 62 
7 11 6.8 7,1 297 279 32 181 a5 53 5 37 106 3.17 «+65 a6 
126 HOLIDAY LAKE APPO"ATTUOX VA 
372339 74382) HOLIDAY LA®2 tite 
lf 4/24 18 s8ha "81 S3 
i« RALNDOMS 4.5 1) Sjee 
*) 17 r.% 8 &e?? 30u 26% a3 143 59 a7 5 7 1 9.39 65 62 
5 le 7.1 8.42 333 3146 5 92 51 4) 5 43 123 1 3 32 b2 
2707 UNSASE Eee APHESST VA 
3732395 PoOdes ABN THST 2158 
H1yY S724 See, 32 Sk? S3 
VI Si R&LINoH: 4,5 lildes 
9 %Y9 %F.3 419 S58 529 bi 188 7 3 5 78 187 0.5 oA 49 
4793 COL SULPA SP BE VES ROCKSBRIDGE VA 
375939 749391) SILLBORO 4408 
sisf 5/15 2.22 jna BRS S53 
ist. SA IND: ‘5.5 lJ lbee 
0 ly 7e 1 7. 5d lb. 127 iu Zi3 48 te “ 11 125 1.9 1s $58 
j lo 7.1 7. 3? 176 te? iu 237 44 2’ ; DD? 122 ée 3 a5 Ha 
/7¥ 











FLD LAB FEP DEP 
D T pH pH ALK ALK COND Al Ca ci cir £& Mg Mo Ne 50, 
2709 MOUNTAIN LAKE 3ILES vd 
372135 89321) EGGLESTON lid2s 
bis 8/25 6.08 22ha 88S S3 
DISTS RAIWoH: 3.5 016ea8 
0 8 6.6 3.33 186 157 27 c) 65 13 3 19 32 0.9 “ a1 
9 6 6.8 6.91 1393 88 21 10 68 17 3 21 69 0.90 7 89 
2719 SUODY PORK C&SEEK TAZEWELL vA 
371520 612015 BSASWELL 7148 
#1/ 4/25 68, ja Bas Sj 
DIST4 RAIHpe: 45 £=%311888 
9 +4109 6.9% 3233 333 294 51 21 130 43 0 30 701 #O.9 21 195 
2777 HUNGRY SOTHESR LAK® S*yvtTi vA 
365215 813125 mABION 671s 
ely 4/25 1.98 SOha BR2 S3 
DIST ®@AINpDH: 4.5 1165008 
09 46 6&6 826 294 26% 53 #911 92 37 3 71 4 9.9 79 117 
a 6 13 6.9 %$%$3453 273 249 51 192 64 ay 3 686 143 0.8 68 106 
7 27912 CLEAR CREZK LASE WASHINGTON VA 
) 363959 629723 WYNDALE 6o1e8 
#1/ 4/25 1.98 20ha 884 53 
DIST} BAINpH: 4.5 116808 
— 17 6.2 8.37 3003 2935 368 25 2011 25 1 33 0 18990 «6242 828 4605 
— e % 115 6.8% 8.36 3209 3092 338 23 2094 a3 1 32 768 2.3 392 813 
.- 
= 2713 CABIN CREEK SRAYSON vA 
P— 16362) BVS129 PARK 10252 
Ww ols 8/25 e8#e a Px2 $3 
DIST} PAtWNoH: 4.5 21988 
nf 4 4.5 6,47 43 s9 21 43 31 15 “ 43 36060. 67 
2715 SeAN LAKE PITTSYLVAS WA 
364215 792014 BLATIRS 198ea 
s1/ 4/26 O.38 Shea 8R1 SS 
DISTS PAIN pd: 4.25 #£=DViistee 
9 1 6.9 7,45 492 353 63 258 77 116 fs] 106 213 2.6 162 137 
9 6.6 6.7% 392 339 71. «259 64 192 6 106 181 8.5 134 163 
2716 COMNZR LAF®e MALIPAL vA 
365529 74461) CONNER LAKZ llje 
417 &/2H OC 28 Sinha BRI S3 
MISTS S8AINDH: 4.5 11} l8#ee8 
9 WW %%9 #%B.93 46) o21 69 196 67 92 3 66 155 46.6 73 252 
3 15 &¥Y $7.83 S51 1019 09 185 73 71 3 73 155 5.6 385 238 
2717 SLAGLES LAK? SREENSVILL VA 
3H9495 773742 BAPORTA jha 
wolf 4&/lo 2698 Saha 8R3 S3 
MTST) PALWOA: 4.6 ll léee 
9 18 72 816 SH 470 99 52 163 45) 618 47 709 06443 +66 333 
3 HW %F.3 3.3% 51D 479 Au 5d 93 le 18 50 120 #O.9 4&3 349 
S19 











FLD LAB FEF DE? 





D T ph Le ALK ALY COND Al Ca Ci CLR ¥ Mg Mn Nea Si 
271d LAKE PRTSC® ISL® OF wT VA 
164539 Jeet05 ete 0508 ée 
417 ast eee, Sha BFe $3 
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APPENDIX B. HISTORICAL DATA 


The historical listing can be separated by state according 
to the folowing code, which appears in the sececr. column of each 
entry, i.e., the second digit of the station code number: 


l. New York 6. West Virginia 
2. New Jersey 7. Virginia 

3. Pennsylvania 8. North Carolina 
4. Delaware ¥. Tennessee 

5. Maryland 


Further information on the specific waters can be found in 
Appendix A. 


Key to column headings for Appendix B: 


1-4 Station numbers (see above) 
6-27 water body name 
29-36 date (year, month, day) of historical sample 
40 depth sample taken 
43-45 pH of sample 
47 method code for pH: l=eiectrometric, 
2=colorimetric 
50-53 alkalinity of sample 
55 method code for alkalinity: l=fixed endpoint, 
2=double endpoint 
58-59 temperature at sample depth (°C) 


A large number of cooperators in state agencies and other sources 
provided the historical data; some of these did not want themselves 
or their data specifically identified. Specific queries concerning 
data for parcicular waters will be accepted and passed cn by the 
authors to the appropriate source. 
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APPENDIX C. STATIONS DROPPED FROM CHEMICAL ANALYSIS DISCUSSION 


The stations listed below were omitted from the data used for 
analysis of chemical relationships for various reasons which are 
coded as follows: 


- poor ionic balance > +100% disagreement 

poor theoretical/observed conductivicy balance 
chloride concentration > 500 yeq/1l 
conductivity > 200 uS/cm 

double endpoint alkalinity > 1000 peq/1l 

color > 30 Pt-Co units 
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Elimination Code 
Station 





a b Cc d P f 





2101 x x x 
2129 x x 
212i 
2123 
2125 
2128 
2132 
2140 
2144 x 

2151 x 
2152 x x x 
2201 x x x 
2204 x 
2206 

2210 

2211 

2214 x 
2215 x 
2224 x 

2334 x 

2337 x 

2343 x 

2349 x 

2350 x 

2352 x 
2353 x x 
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x KK XK 
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(continued) 





Elimination Code 








Station 
a b Cc d e f 
2354 x 
2362 x x x 
2363 x x 
2365 x 
2374 x 
2376 x 
2377 x 
2406 x x 
2411 x 
2503 x 
2507 x 
2510 x 
2511 x 
m 2513 x x 
2522 x 
2523 x 
2524 x 
2525 x 
2603 x x 
2627 x 
2628 x 
2630 x 
2634 x x 
2635 x 
2637 x 
2640 x 
2701 x 
2702 x x 
2712 x 
2718 x 
2803 x 
2903 F x 
2905 x 


2906 x x 








The net 


£ stations 


218; 


a 22% decrease. 


result of 


used 


New York ll of 
New Jersey 8 of 
Pennsylvania l4 0 
Delaware 2 of 
Maryland 9 of 
West Virginia 8 of 
Virginia 4 of 
North Carolina l of 
Tennessee 3 of 
Total 60 of 


these 
for discussion of 
Removals by state 


removals 
the 
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is to decrease the number 
chemical data from 278 


were: 

46 (24%) 
31 (26%) 
88 (16%) 
21 (10%) 
24 (38%) 
23 (35%) 
24 (17%) 
1® (10%) 
ll (27%) 
278 (22%) 
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| the waters were acidified. Few of the waters showed high concentrations of 

aluminum, and there was only one case of a potentially toxic aluminum/pH 
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As the Nation's principal conservation agency. the Department of 
the interior has responsibility for most of our nationally owned 
public lands and natural resources. This includes fostering the 
wisest use of our land and water resources. protecting our fish and 
wildiife. preserving the environmental and cultural values of our 
national parks and historical places, and providing for the enjoy- 
ment of life through outdocr recreation The Department assesses 
our energy and mineral resources and works to assure that their 








development is in the bes: interests of ali our people. The Depart- f Ls, \ 
FISH & WILDLIFE 
ment also has a major responsibility for American indian reservation SERVICE 





communities and for people who live in island territories under US 
administration 
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